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Abstract 


This  report  summarizes  the  work  performed  during  the  second  year 
of  an  AFOSR  sponsored  research  program  that  was  primarily  concerned  with 
the  development  of  an  analytical  technique  for  determining  the  radiated 
sound  field  from  axisymmetric  jet  engine  inlet  configurations.  The 

analytical  technique  employed  is  based  upon  an  integral  representation  of 
the  external  (radiation)  solutions  of  the  Helmholtz  equation  which  describe 
the  sound  fields  external  to  a given  body  under  either  no  flow  or  constant 
velocity  flow  situations.  The  integral  representation  developed  during  the 
course  of  this  research  program  is  different  from  earlier  works  in  the  sense 
that  it  not  only  yields  the  correct  (unique)  solution  for  ail  radiation  problems 
at  all  frequencies,  but  that  the  resulting  integral  equations  contain  no  strong 
(i.e.,  non-integrable)  singularities  and  therefore  can  be  solved  by  straight 
forward  numerical  techniques.  As  part  of  this  research  effort  two 

extremely  flexible  computer  programs  were  developed  for  the  solution  of 
these  axisymmetric  integral  equations.  These  programs  can  be  used  for  any 
closed  axisymmetric  body  with  any  combination  of  boundary  conditions  on 
its  surface,  without  any  modification  of  the  computer  codes,  to  provide 
accurate  solutions  for  the  acoustic  properties  (i.e.,  the  acoustic  pressure, 
normal  velocity,  and  admittance)  both  on  the  body  itself  and  anywhere  in  the 
field  surrounding  the  body.  Also,  as  part  of  this  research  effort,  experiments 
are  now  being  conducted  with  various  axisymmetric  configurations  t ' 
provide  the  data  that  will  be  used  to  check  the  validity  of  the  theoretical 
predictions.  Some  preliminary  experimental 
report. 
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I.  Introduction 

This  report  summarizes  the  results  obtained  during  the  second  year 
of  support  under  AFOSR  contract  number  F49620-77-C-0066.  This  contract 
was  initiated  on  February  1,  1977  and  the  results  obtained  during  the  first 
year  of  support  are  contained  in  AFOSR  Interim  Scientific  Report  AFOSR- 
TR-78-0696. 

The  main  objective  of  the  research  program  conducted  under  this 

contract  was  to  develop  an  analytical  technique  for  predicting  the  sound 

field  radiated  from  axisymmetric  jet  engine  inlet  configurations  and  to 

check  the  validity  of  these  predictions  with  relevant  experimental  data.  The 

development  of  this  analytical  solution  technique  was  motivated  by  the  need 

for  a theoretical  approach  that  could  be  used  to  predict  the  effects  of  sound 

source  modifications  and  of  sound  suppression  devices  (such  as  acoustically 

lined  surfaces  and  splitters)  upon  the  sound  field  radiated  from  the  inlet 

without  having  to  resort  to  costly,  full  scale  experimental  testing.  The 

experimental  investigations  are  necessary  not  only  for  comparison  with  the 

results  of  the  analytical  technique  (which  has  shown  extremely  good 

.12  3* 

agreement  with  known  exact  solutions)  ’ ’ but  also  to  assist  in  the 
determination  of  the  correct  analytical  form  for  describing  the  boundary 
conditions  necessary  to  accurately  model  sound  suppression  materials  in  the 
computer  programs. 


These  references  were  included  as  appendices  of  the  aforementioned 
Interim  Scientific  Report,  AFOSR-TR-78-0696. 
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During  the  first  year  of  this  contract  the  integral  equations  and 

computer  programs  required  for  investigating  the  sound  fields  radiated  from 

certain  simple  geometries  were  developed.  Specifically,  two  different 

formulations  of  the  problem  were  developed;  that  is,  the  full  three 

dimensional  formulation  and  the  axisymmetric  formulation.  Employing  both 

these  formulations  the  sound  radiation  from  both  a sphere  and  a finite  length 

cylinder  were  investigated.  In  these  studies  the  effects  of  different 

boundary  conditions  on  the  accuracy  of  the  integral  solution  technique  was 

investigated  by  comparing  the  predictions  of  this  approach  for  the  sound 

fields  radiated  by  simple  sources  with  available  exact  solutions  obtained 

using  the  Seperation  of  Variables  Technique  (which  is  only  applicable  to 

simple  geometries  such  as  those  investigated).  This  work  is  described  in 

detail  in  three  publications  ’ ' (See  footnote  on  the  previous  page.)  and  it 

ARC** 

was  presented  at  three  conferences  ’ ’ . In  summary  the  first  year  of 

study  has  shown  that  the  integral  solution  technique  developed  under  this 
contract  is  both  accurate  and  efficient  from  a numerical  point  of  view. 

During  the  second  year  of  this  contract  two  general  computer 
programs  were  developed  for  the  determination  of  the  acoustic  fields  both 
on  and  around  general  axisymmetric  bodies  with  general  boundary 
conditions.  It  will  be  noted  here  that  although  the  geometry  of  the  bodies 
under  consideration  are  constrained  to  be  axisymmetric  the  allowable 
acoustic  modes  are  not;  that  is,  any  cylindrically  symmetric  acoustic  mode 
which  may  be  present  can  be  solved  for  using  these  computer  programs. 
These  programs  were  written  in  Extended  Fortran  IV  and  they  are  presented 


Letters  refer  to  the  various  conference  presentations  conducted  in 
connection  with  this  program.  These  presentations  are  listed  in 
Appendix  F. 
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in  Appendix  A.  They  have  been  fully  checked  out  on  the  Georgia  Tech  CDC 
Cyber  70/74  computer  and  have  been  used  to  theoretically  predict  the 
radiated  sound  fields  that  are  associated  with  both  lined  and  unlined  straight 

tj, 

ducts  and  an  actual  jet  engine  inlet  configuration  shown  in  Fig.  1. 

The  details  and  results  of  these  investigations  are  described  in 
References  5 and  6 which  are  reproduced  in  Appendices  B and  C 
respectively.  The  latter  of  the  two  papers  has  also  been  accepted  for 
presentation  at  the  AIAA  5th  Aeroacoustics  Conference  in  March  1979  in 
Seattle,  Washington.  Probably  the  most  significant  results  of  this 
investigation  is  the  result  that  the  optimum  theoretical  admittance  values 
found  for  reducing  the  sound  radiated  from  a straight  duct  do  not  necessarily 
have  the  same  effect  in  a dimensionally  similar  (i.e.,  the  same  length  to 
diameter  ratio)  inlet.  That  is,  the  geometrical  details  of  an  axisymmetric 
body  need  to  be  taken  into  consideration  when  calculating  optimum 
admittance  values  for  acoustic  liners.  Another  resuit  of  this  investigation  is 
that  the  admittance  at  the  entrance  plane  of  a straight  duct  or  an  inlet  is 
not  constant  in  the  radial  direction.  The  assumption  of  a constant 
admittance  value  at  the  entrance  plane  of  a duct  is  common  to  many  current 
theoretical  analyses  of  the  duct  radiation  problem  and  thus  can  be  a large 
source  of  error  in  these  analyses. 

Since  the  use  of  a constant  admittance  value  at  the  entrance  plane 
of  a duct  is  so  common  a series  of  computer  runs  were  performed  to 
determine  just  how  accurate  this  assumption  is.  The  runs  were  done  for  a 
straight  duct  configuration  and  various  admittance  values  were  calculated 
and  compared  with  certain  "classical  values".  The  details  of  this  analysis  and 
the  results  are  presented  in  Appendix  E.  In  Appendix  E a draft  copy  of  a 
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technical  note,  to  be  submitted  for  publication  in  the  A1AA  Journal,  is 
reproduced. 

Part  of  the  efforts  expended  under  this  contract  during  the  second 
year  of  study  were  devoted  to  the  development  of  experimental  data  that 
could  be  used  to  check  the  validity  of  the  accompanying  analytical  studies. 
These  efforts  consisted  of  the  design  and  development  of  appropriate 
experimental  setups  and  the  conduct  of  the  required  experiments,  which  are 
currently  in  progress.  The  experimental  data  acquired  under  this  program  is 
being  compared  with  corresponding  theoretical  predictions.  Some 
preliminary  results  of  these  comparisons  are  now  available  and  they  are 
presented  in  Section  IV  of  this  report.  The  work  conducted  during  the 
second  year  of  this  contract  which  is  not  presented  in  the  papers  contained 
in  Appendices  B and  C is  summarized  in  the  following  sections. 


I 
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II.  Analytical  Effort 


The  analytical  efforts  during  the  second  year  of  this  AFOSR 
contract  consisted  of  the  development  of  the  general  computer  programs  for 
handling  axisymmetric  geometires  and  their  use  to  perform  certain 
parametric  studies  of  interest.  The  details  and  results  of  these  studies  are 
presented  in  Appendices  B,  C,  and  E,  and  the  computer  programs  themselves 
are  presented  in  Appendix  A.  The  programs  were  developed  so  that  they 
could  not  only  be  used  for  parametric  analytical  studies  but  also  so  that 
they  could  easily  accept  experimental  data  for  any  configuation  without 
any  major  changes  in  the  computer  code.  Another  part  of  the  analytical 
effort  was  concerned  with  the  determination  of  the  admittance  of  the  liner 
which  was  used  in  the  experimental  phase  of  this  program  and  with  the 
redesign  of  the  liner  for  future  testing  (See  Appendix  D.). 

* The  initial  testing  under  this  program  was  conducted  with  an 

available  acoustic  liner  which  was  developed  in  a related  combustion 
instability  program.  This  liner  was  tuned  for  maximum  damping  at  a 
frequency,  ~ 7W  Hertz,  which  is  above  the  first  tangential  (i.e.,  IT)  mode, 
~695  Hz  ,of  the  duct  under  investigation.  Since  we  are  mainly  interested  in 
runing  experiments  below  the  IT  cut-off  frequency  to  facilitate  both  the 
data  reduction  and  the  comparison  of  experimental  and  analytical  values  all 
of  the  preliminary  testing  was  done  in  the  frequency  range  300-650  Hz 
which  is  below  the  IT  cut-off  frequency.  Thus,  the  liner  was  not  expected 
to  exert  much  attenuation.  Since  most  of  the  planned  future  testing  will  also 
be  conducted  over  a frequency  range  below  the  IT  mode  of  the  duct,  the 
available  liner  will  be  retuned  to  be  more  effective  below  the  IT  cut-off 

ft 

frequency.  The  retuning  of  the  liner  is  discussed  in  detail  in  Appendix  D. 
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III.  Experimental  Investigations 


The  main  objective  of  the  experimental  phase  of  this  program  was  to 
obtain  experimental  data  that  could  be  readily  compared  with  the 
predictions  of  the  analytical  models  developed  under  this  program. 
Specifically,  the  sound  fields  radiated  from  lined  and  unlined  axisymmetric 
duct  configurations  were  to  be  measured  and  compared  with  corresponding 
theoretical  predictions.  Since  these  studies  were  concerned  with  the 
measurements  of  the  radiated  sound  fields,  all  of  the  required  experiments 
were  conducted  in  an  anechoic  chamber  whose  properties  are  described  in 
Fig.  2.  A typical  experimental  setup  utilized  in  the  course  of  this  study  is 
shown  in  Fig.  3;  it  consists  of  a lined  axisymmetric  duct  with  a sound  source 
at  one  end  and  an  open  termination  at  the  other  end.  The  test  body  (i.e.,  the 
inlet)  is  placed  on  one  side  of  the  anechoic  chamber  and  an  array  of 
microphones  is  used  to  measure  the  radiated  sound  field.  The  latter 
consisted  of  5 Bruel  and  Kjaer  condenser-type  microphones.  The  sound 
source  was  a University  driver  and  it  was  placed  at  the  throat  of  a nozzle 
that  was  connected  to  the  axisymmetric  body  tested  as  shown  in  Fig.  **.  The 
available  liner  used  for  the  preliminary  testing  was  not  specifically  designed 
for  maximum  effectiveness  in  the  frequency  range  where  most  of  the  testing 
was  done.  It  was  tuned  for  maximum  damping  at  a frequency  higher  than  the 
IT  mode  of  the  duct  under  study.  It  consisted  of  180  Helmholtz  resonators 
(20  radial  rows  by  9 axial  rows)  which  were  closed  off  during  the  hard  walled 
testing  (See  Appendix  D for  the  admittance  calculation).  A diagram  of  the 
driver-nozzle-liner  set-up  and  of  one  of  the  Helmholtz  resonators  appears  in 
Fig.  4. 
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ANECHOIC  CHAMBER 

PLAN  VIEW 


Figure  3.  Experimental  Set-up  for  Radiation  Studies. 
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This  experimental  set-up  was  found  to  work  very  well  in  that  the 
acoustic  waves  generated  by  the  driver  remained  essentially  plane  up  to  the 
nozzle-liner  interface.  This  was  determined  by  traversing  a 1/4"  B<5cK 
microphone  radially  across  the  nozzle-liner  interface.  This  being  the  case 
only  a single  measurement  of  the  acoustic  pressure,  which  is  needed  for 
input  to  the  analytical  computations,  was  taken  at  this  plane  for  each  test 
condition. 
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IV.  Preliminary  Experimental  Results  and  Comparisons 


Some  preliminary  experiments  have  been  conducted  with  the 
straight  duct  configuration  shown  in  Figs.  3 ad  4.  Both  lined  and  unlined 
configurations  were  tested;  however,  not  much  difference  was  noted  as  all 
the  tests  were  conducted  at  frequencies  below  the  first  tangential  (i.e.  IT) 
mode  of  the  duct  and  therefore  in  a region  where  the  liner  is  only  marginally 
effective  as  discussed  in  the  previous  sections  (Also  see  Appendix  D.).  Tests 
were  conducted  in  the  frequency  range  300  to  650  Hertz  with  50  Hz 
increments.  The  microphones  in  the  field  were  placed  on  a circular  arc  with 
a radius  of  seven  feet  centered  at  the  entrance  plane  of  the  duct  (See  Figure 
3.).  The  microphones  were  placed  at  increments  of  11K°  from  the 
centerline  of  the  duct  to  90° . 

Comparisons  between  the  experimental  results  and  the  theoretical 
results  were  made.  Since  an  experimentally  measured  value  of  the  sound 
pressure  level  at  the  nozzle-duct  interface  was  used  as  input  for  the 
computer  programs,  the  accuracy  of  the  calculated  far  field  can  only  be 
expected  to  be  as  good  as  this  measurement.  Other  sources  of  error  are  the 
"imperfections"  of  the  anechoic  chamber  at  various  frequencies  as  shown  in 
Fig.  2,  instrumentation  errors,  and  the  lack  of  perfect  correspondence 
between  the  experimental  and  analytical  configurations.  In  this  connection  it 
should  be  pointed  out  that  the  theoretical  model  employs  a spherical 
termination  at  the  rear  end  of  the  straight  duct  as  shown  in  Fig.  5.  This  was 
done  to  improve  the  efficiency  of  the  computer  programs  as  it  has  been 
shown  through  theoretical  studies  that  the  shape  of  the  termination  of  the 


duct  exerts  little  influence  upon  the  acoustic  field  in  the  forward  half  plane 
of  the  duct. 

Comparisons  between  the  measured  and  calculated  results  shows 
good  qualitative  agreement  for  both  the  hard  walled  (i.e.,  See  Fig.  6.)  and 
the  lined  (i.e.,  See  Fig.  7.)  duct  configurations.  Good  agreement  is  observed 
between  the  hard  walled  and  lined  wall  cases  in  that  the  errors  follow  the 
same  patterns  (i.e.,  Compare  Figs.  6 and  7.)  which  suggest  that  most  of  the 
observed  errors  are  due  to  the  "non-anechoicness"  of  the  anechoic  chamber. 
The  measured  data  will  be  further  analyzed  in  the  future  and  they  will  be 
published  together  with  additional  data  collected  during  the  next  year  of 
study  under  this  contract. 
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Appendix  A 


Computer  Programs  Developed  for  this  Project 

The  computer  programs  described  in  this  appendix  were  developed 
during  the  second  year  of  this  AFOSR  research  project.  They  are  very 
general  in  that  they  can  be  used  to  find  the  acoustic  properties  (i.e.,  the 
acoustic  potential,  normal  acoustic  velocity,  and  the  admittance)  both  on 
the  surface  and  surrounding  any  finite  (closed)  axisymmetric  body.  Also,  the 
boundary  conditions  can  be  specified  at  each  point  on  the  surface  of  the 
body;  the  only  restriction  is  that  the  admittance  may  not  be  specified 
everywhere  on  the  body  as  then  the  solution  of  the  governing  integal 
equation  is  non-unique. 

The  first  computer  program  developed  for  this  project  calculates  the 
surface  distributions  of  the  acoustic  quantities  for  a general  axisymmetric 
configuration.  Required  inputs  for  this  program  are  the  geometric 
description  of  the  body,  the  boundary  conditions,  and  the  problem 
specification  data  (i.e.,  the  mode  number  of  the  acoustic  wave  being  solved 
for  and  the  wave  number).  The  geometric  input  data  required  includes  the 
p - z coordinates  of  the  points  in  the  center  of  each  integration  interval,  P, 
the  two  integration  points  on  either  side  of  the  center  point,  Q,  the  normal 
to  the  body  at  the  center  point,  NORMAL,  and  the  length  of  each 
integration  interval  in  the  p -z  plane,  LENGTH.  Since  the  body  is  assumed  to 
be  axisymmetric  the  program  takes  care  of  choosing  the  integration  points 
in  the  0 direction  through  the  use  of  a 96  point  Gaussian  integration 
formula. 
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The  boundary  conditions  are  specified  in  two  vectors.  The  first 
vector,  ICHECK,  specifies  what  type  of  boundary  condition  is  known  at  each 
point  in  the  p-z  plane. 


ICHECK 


<p  known 

d cp/dn  known 
y known 


(A-l) 


The  second  vector,  C1NDATA,  contains  the  actual  complex  value  of  the 
boundary  condition  at  each  point.  Finally  the  problem  specification  data 
consists  of  specifying  the  wave  number  k and  the  mode  number  m (i.e.  m = 0 
is  the  axisymmetric  mode.) 

This  program  prints  out  all  the  input  data  and  all  the  acoustic 
quantities  on  the  surface  of  the  body.  It  also  creates  an  output  file  which  is 
read  by  the  next  program  to  calculate  the  acoustic  qualtities  in  the  field 
surrounding  the  body. 


■ 


c 

c 


PROGRAM  EXPCYL  ( INPUT,  OUTPUT,  TAPE 10,  TAPE 1 1 , 
TAPE5  = INPUT,  TAPE6  = OUTPUT) 


C*  * 

C*  * 

C*  THIS  PROGRAM  CALCULATES  EITHER  THE  ACOUSTIC  POTENTIAL  OR  THE  * 

C*  ACOUSTIC  VELOCITY  OR  THE  EFFECTIVE  ADMITTANCE  ON  THE  SURFACE  * 

C*  OF  ANY  AXISYMMETRIC  BODY  USEING  THE  METHOD  OF  BURTON  S MILLER  * 

C*  WITH  MY  INTERPRETATION  OF  THE  MOST  SINGULAR  COMPONENT.  * 

C*  * 

C****«****************************************************************** 

c*  * 

C*  WILLIAM  A.  BELL’S  OPTIMAL  VALUE  OF  ALPHA  ( = I/K)  IS  EMPLOYED.  * 

C*  * 

c*********************************************************************** 

c*  * 

C*  A CYLINDRICAL  FORMULATION  OF  THE  PROBLEM  IS  EMPLOYED.  * 

C*  * 

c*  * 

c*******************^************************************************** 


c 

c 


c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


REAL  K,  LENGTH,  NORMAL 

COMPLEX  ALPHA,  CEQN,  CEXACT,  CPHI,  CVEL.  CY,  IK,  IKSO,  TWOPIA, 

, C INDATA 

COMMON  /I/  M,  NP.  NPP 1 
COMMON  /XV  K,  PI,  TWOPI 
COMMON  /C/  ALPHA.  IK,  IKSO,  TWOPIA 
COMMON  SID/  I CHECK  (102) 

COMMON  /RD/  LENGTH  (102).  NORMAL  (102,  2),  P (102,  2). 

O ( 102,  2.  2) 

COMMON  /CD/  CEON  (102,  103),  CEXACT  (102),  CPHI  (102),  CVEL  (102), 
CY  (102),  CINDATA  (102) 

COMMON  /NGAUSS/  NCAUST 
COMMON  /GAUSS/  CAUSNT  (48,  2) 

CALL  INPUT 

INTEGER  CONSTANTS 

NPP 1 = NP  + 1 

REAL  CONSTANTS 

TWOPI  = 2.0  * PI 

COMPLEX  CONSTANTS 

IK  = (0.0,  1.0)  * K 

IKSO  = IK  * IK 
ALPHA  = (0.0.  1.0)  / K 

TWOPIA  = ALPHA  * TWOPI 


C 

C 

C 

C 


CALL  F.QN 
CALL  GAUSS 
CALL  OUTPUT 
STOP  "NORMAL* 
END 


20 


- * — r 

I MMM  " - — 


BLOCK  DATA 
REAL  K 


COMMON  /I/  M,  NP,  NPP1 
COMMON  /R-  K,  PI,  TWOPI 
COMMON  /NCAUSS/  NGAUST 
COMMON  /GAUSS/  GAUSNT  (48,  2) 


DATA  NP  / 102  / 

DATA  PI  / 3.1415926535898  / 
DATA  NGAUST  / 40  / 


DATA 


DATA 


DATA 


DATA 


( ( GAUSRT( I , J),  J = 
/ 0.01627674484960, 
0.04801298513605, 
G.00129749546443, 
0. 1 136950501 1067, 
0. 14597371465490, 
0. 17809600236762, 
0.21003131046057, 
0.24174315616384, 
0.27319801259105, 
0 . 30436494435450 , 
0 . 33520852209263 , 
0.36569686147231, 
( ( GAUSNT(  I , J) , J = 
/ 0.39579764982891, 
0.42547890040730, 
0.45470942216774, 
0 . 48345797392060 , 
0.51  169417715467, 
0.53938010032436, 
0.56651041856140, 
0.59303236477757. 
0.61092504012547, 
0.64416340370497, 
0.66871031004392, 
0.69256453664217, 
( ( GAUSNTC I , J) , J = 
/ 0.71567681234097, 
0 . 73003064374440 , 
0.759602341 17665. 
O . 78036904306743 , 
0.80030874413914, 
0.01940031073793, 
0.03762351 1220 19, 
0 . 05495903343460 , 
0.07130850590930, 
0.88689451740242. 
0.90146063531585, 
0.91507142312090, 
( ( GAUSNTC I . J ) . J = 
/ 0.92771245672231, 
0.93937033975276. 
0.95003271778444, 
0.95968829144074, 
0 . 96832682846326 . 
0.97593917458514, 
0.98251726356301  , 
0.98805412632962, 
0 . 99254390032376 . 
0.99598184298721 , 
0.99036437586318. 
0.99960950388323, 


I,  2),  I * 1,  12) 
0.03255061449236, 
0.03251611071387, 
0.03244716371406, 

0 . 03234382256858 , 

0 . 03220620479403 , 
0.03203445623199, 
0.03182875889441 , 
0.03158933077073, 
0.03131642559686, 
0.03101033258631 , 
0.03067137612367, 

0. 03029991542083  / 

1,  2).  I = 13,  24) 
0.02989634413633, 
0.02946108995817, 
0.02899461415056, 
0.02849741 106509, 
0.0279700076  1685. 
0.02741296272603, 
0.02682686672559, 
0.02621234073567. 
0.02557003600535, 
0.02490063322248, 
0.02420484179236, 

0. 02348339908593  / 

1,  2),  1 = 25,  36) 
0.02273706965833, 
0.02196664443874, 
0.02117293989219, 
0.02035679715433, 
0.01951908114015, 
0.01866067962741, 
0.01778250231605, 
0.01688547986425, 
0.01597036290256, 
0.01503872102699, 
0.01409094177231 , 

0. 01312822956696  / 

1.  2),  I - 37,  48) 
0.01215160467109, 
0.011 16210209984, 
0.01016077053501 , 
0.00914067123078, 
0.00812687692570, 
0.007096470791 15, 

O . 00605854550424 , 
0.00301420274293, 
0.00396455433844, 

0 . 0029 1073181793, 
0.00185396078895, 
0.00079679206555  / 


c 


SUBROUTINE  »NPUT 


C 


c 

c 

c 

c 

c 

100 

c 


REAL  K,  LENGTH.  NORMAL 

COMPLEX  CEQN,  CEXACT.  CPHI.  CVEL,  CY.  C INDATA 

COMMON  /\/  M.  NP,  NrPl 
COMMON  /XV  K.  PI,  TWOP [ 

COMMON  /'ll)/  I CHECK  (102) 

COMMON  /IU)/  LENGTH  (102),  NORMAL  (102,  2),  P (102,  2), 

Q ( 102,  2,  2) 

COMMON  /CD/  CEQN  (102,  103).  CEXACT  (102),  CPHI  (102),  CVEL  (102) 

CY  (102),  CINDATA  (102) 


READ 

(5. 

100) 

( ( P( I , J)  , 

J = 1,  2), 

I = 

1,  NP) 

READ 

(5, 

100) 

( ( ( Q(  I , J , 

KK)  , KK  = 1 

. 2) 

, J = 1,  2) 

READ 

(5. 

100) 

( ( NORMAL(  I 

, J)  , J = 1, 

2)  , 

I = 1,  NP) 

READ 

(5, 

100) 

( LENGTH!  I) 

, I = 1,  NP) 

FORMAT  (80 10.0) 


PRINT  101 


C 

101 


C 


C 

102 


C 

C 

104 

C 


FORMAT  ( " 1GE0MF.TR I C INPUT  DATA:"  /// 

" ",  25X,  "-Q( RHO,  Z) " / 

« ",  4X,  "N'\  2 IX,  "P(RHO,  Z)’.  38X,  "NORMAL(  RHO,  Z)*. 

27X,  "LENGTH"  / 

» ",  25X,  "+Q(RHO,  Z) " //  “ ") 


WRITE  (6,  102)  (Q(I,  1,  1),  Q(I,  I,  2),  I,  P( I , 1),  P( I , 2), 

NORMAL! 1 , 1) , NORMAL! I,  2) , LENGTH!  I) , Q(  I , 2, 
Q( I,  2,  2),  I = I,  NP) 


FORMAT  ( " 

N 


II 


15X,  "( “,  F 13. 10, 

2X,  13,  10X,  "(“. 

» ( II 

I5X,  ”(  ”,  F13. 10, 


",  ”,  F13.  10,  ">  " ✓ 

F 13 . 10,  ”,  ",  F13. 10. 
F13. 10,  ”,  ",  F13. 10, 

F13.  10  ✓ 

",  ",  F 13.  10,  ”)  " / " 


") 


2 IX, 
17X, 


READ  (5,  104)  (ICHECK(I),  I = 1 , NP ) 


FORMAT  (1615) 


READ  (II.  1 100)  M,  K 


C 

1100  FORMAT  (120,  C20.0) 


C 

READ  (11,  100)  (CINDATA(I),  I = I . NP) 

C 

RETURN 

END 


1) 


SUBROUTINE  EON 


i 

* 


C 


c 


c 

c 

c 


2 


c 


3 


4 


6 

1 

C 

C 

C 


C 

9 

U 

7 

C 


REAL  K,  LENGTH , NORMAL,  NP1,  NP2,  NRHOP , NRHOQ,  NZP,  NZQ 
COMPLEX  ALPHA,  CEQN.  CEXACT,  CPHI,  OTX.  CY,  IK,  IKBQ, 

TWOPIA,  C1NDATA 

COMMON  /!/  M.  NP,  NPPl 
COMMON  /R/  K,  PI,  TWOP I 
COMMON  /C/  ALPHA,  IK,  IKSQ.  TWOPIA 
COMMON  /ID/  ICHECK  ( 102) 

COMMON  /RD/  LENGTH  (102),  NORMAL  (102,  2),  P (102,  2). 

Q ( 102,  2,  2) 

COMMON  /CD/  CEQN  (102,  103),  CEXACT  (102),  CPHI  (102),  CVEL  ( 
CY  (102',  CINDATA  (102) 

1 N ITALIZE  MATRIX 

1)0  I I = 1,  NP 

DO  2 .1  = 1 . NPPl 

CEQN!  I,  J)  « (0.0,  0.0) 

CONTINUE 
PI  •=  P(  I . I) 

P2  = P(  1 , 2) 

NP 1 a NORMAL! I.  I) 

NP2  = NORMAL!  1 . 2) 

IP  (I  CHECK!  I ) ) 3 . 4 , 5 
CONTINUE 

CEXACT!  1)  = CINDATA!  I ) 

CVEL! I ) a (0.0,  0.0) 

CY! I)  = (0.0,  0.0) 

CEQN! I , NPPl)  = TWOP 1 * CEXACT! I ) 

CEQN!  I . 1)  = -TWOPIA 

CO  TO  6 

CONTINUE 

CEXACT!  I ) = (0.0,  0.0) 

CVEL! I ) = CINDATA! 1) 

CY! I)  = (0.0.  0.0) 

CEQN!  I,  NPPl)  a TWOPIA  * CVEL( I) 

CEQN! I,  I)  = -TWOP I 

GO  TO  6 

CONTINUE 

CEXACT(I)  = (0.0,  O.O) 

CVEL! I ) = (0.0,  0.0) 

CY!  I ) = CINDATA!  I ) 

CEQN! I,  I)  = -TWOP I * (1.0  + ALPHA  * CY(I)) 

CONTINUE 

CONTINUE 

XI  INTEGRATION 

DO  7 J = 1,  NP 
NRHOQ  = NORMAL!  J,  1) 

NZQ  = NORMAL!  J , 2) 

GAUSZ  = LENGTH!  J)  * PI 
DO  O I = 1,  NP 
RHOP  = P(I,  1) 

ZP  a P(  1 , 2) 

NRHOP  = NORMAL!  1 . 1 ) 

NZP  = NORMAL!  I,  2) 

DO  9 L = 1,  2 

RHOQ  = Q(.I,  L,  1) 

ZQ  a 0!.I,  L,  2) 

CALL  CALC  (RHOP,  ZP,  NRHOP,  NZP.  RHOQ.  ZQ,  NRHOQ,  NZQ,  GAUSZ 
I,  .1) 


CONTINUE 

CONTINUE 

CONTINUE 

RETURN 

END 


102)  , 
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SUBROUT IRE  CALC  (RHOP,  ZP,  NRHOP,  NZP,  RHOQ.  ZQ,  NRHOQ,  NZQ, 

GAUSZ,  I,  J) 

REAL  K,  NDOTN,  NRHOP.  NRHOQ,  NRHOQP , NZP,  NZQ,  NZQP 
COMPLEX  ALPHA,  CEQN.  CEXACT,  CPHI.  CVEL,  CY,  FI,  F2,  G,  GP, 

CPP,  IK,  IKSQ,  II,  12,  Kl,  K2,  TWOPIA,  CINDATA 

COMMON  /I/  M,  NP,  NPPI 

COMMON  /O'  K,  PI,  TWOPI 

COMMON  / C / ALPHA,  IK,  IKSQ,  TWOPIA 

COMMON  /ID/  ICHECK  (102) 

COMMON  /Cl)/  CEQN  (102,  103),  CEXACT  (102),  CPHI  (102),  CVEL  (102), 

. CY  (102),  CINDATA  (102) 

COMMON  /NGAUSS/  NGAUST 
COMMON  /GAUSS/  GAUSNT  (48,  2> 


ZD  = ZQ  - ZP 

ZSQ  = ZD  * ZD 

RHOSQ  = (IUIOQ  - RHOP) *#2 

RII0QP2  = 2.0  * RIIOO  * RHOP 

NIUIOQP  = NRHOQ  * NRHOP 

NZQP  = NZQ  * NZP 

THETA  INTEGRATION 

DO  1 IT  = 1,  NGAUST 
THETA  = PI  * GAUSNT(  IT,  I) 

IF  (1  .EQ.  J)  THETA  = PI  * (1.0  - GAUSNT(  IT,  1>> 

GAUSZT  = GAUSZ  * CAUSNT(  IT,  2) 

COST  = COS  (THETA) 

COS Ml  = COS  (M  * THETA) 

R = SORT  ( RHOSQ  + ZSQ  + RH0QP2  * (1.0  - COST) ) 

DRDNO  = ( NRHOQ  * ( RHOQ  - RHOP  * COST)  + NZQ  * ZD)  / R 
DRDNP  = ( NRHOP  * ( RHOP  - RHOQ  * COST)  - NZP  * ZD)  / R 
NDOTN  = NRHOQP  * COST  + NZQP 

G = IUIOQ  * CEXP  (IK  * R)  / R 
GP  * C * ( IK  - ( 1.0  / R)> 

GPP  = G * ( IKSQ  - (3.0  * IK  / R)  + (3.0  / (R  * R)>) 

1 1 -•  c * cosvrr 

12  = ALPHA  * GP  * DRDNP  * COSPTT 
FI  = ALPHA  * G * IKSQ  * NDOTN 

F2  = ALPHA  * (CPP  * DRDNP  * DRDNQ  - ( GP  * NDOTN  / R) ) 

Kl  - GP  * DRDNQ  * COSMT 
K2  = F2  * COSMT 

IF  ( ,CHECK(  I ) -NE.  -1)  CO  TO  3 

CEQN(  I , NPPI)  = CEQN(  I , NPPI)  + GAUSZT  * CEXACT(l)  * (FI  + F2) 


GO  TO  4 
CONTINUE 


CEQN(  I , I)  = CEQN( I , I)  - GAUSZT  * (FI  + F2) 


CONTINUE 

IF  ( ICHECK( J) ) 11,  12.  13 
CONTINUE 

CEURII,  J>  = CEQN(  I,  J)  - CAUSZT  * (II  + 12) 

CF.QIK  I . NPPI)  = CEQN(  I , NPPI)  - GAUSZT  * (Kl  + K2) 

* CEXACT!  J) 


GO  TO  14 
CONTI  NUF. 

CEQN(  I , J)  * CEQN ( I . J)  + GAUSZT  * (K!  + K2) 
CF.QN'l.  NPPI)  = CEQN(  I , NPPI)  + GAUSZT  * ( 1 1 ♦ 12) 

* CVEL(J) 


GO  TO  14 
CONTINUE 

CEQNl  I,  J)  « CEQN(I,  J)  + GAUSZT  * ((Kl  + K2>  - 

CY( J)  * (II  ♦ 12) ) 


CONTINUE 
CONT I nue 


RETURN 
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SUBROUTINE  GAUSS 


IMPLICIT  COMPLEX  ( C) 

C 

COMMON  /I/  MM,  NP,  NPP1 

COMMON  /CD/  CEQN  (102,  103),  CEXACT  (102),  CPHI  (102),  CVEL  (102), 
CY  ( 102) . CINDATA  ( 102) 

C 

C UPPER  TR1ANCUL ARIZE  MATRIX 

C 

DO  1 J = I,  NP 
JP1  = .1  + 1 
CSAVE  = CEQN  (J,  J) 

CEQN( J,  J)  = (1.0.  0.O) 

DO  2 L = JP1,  NPP 1 

CEQIKJ,  L)  = CEQN(J.  L)  / CSAVE 

2 CONTINUE 

IF  (J  .EQ.  NP)  GO  TO  3 
DO  4 M = JP1 , NP 
CSAVE  = CEQN(M,  J) 

DO  5 I = ,IP1,  NPP  1 

CF.QN(  M,  I)  = CEQN(M,  I)  - CEQN(J,  I)  * CSAVE 

5 CONTINUE 

4 CONTINUE 

1 CONTINUE 

C 

C RACK  SUBSTITUTION 

C 

3 CONTINUE 

CSUM  = (0.0.  0.0) 

DO  6 I = 1.  NP 
NPMI  = NP  - I 
NPP 1 Ml  = UPMI  + 1 

CPHI(NPPIMI)  = CEQN(  NPP1MI , NPP1)  - CSUM 
IF  ( I .EQ.  NP)  CO  TO  7 
CSUM  = (0.0,  0.0) 

DO  0 J = I . I 
NPP I MJ  = NPP i - J 

CSUM  = CPU I ( NPP1MJ)  * CEQN( NPMI , NPP1MJ)  ♦ CSUM 
O CONTINUE 

6 CONTINUE 

7 CONTINUE 
C 

RETURN 

END 


— 


c 


c 


c 


c 


c 

c 

lee 


c 


3 

c 

c 

■ ei 

c 


4 

c 

c 

102 

C 


c 


c 

103 

C 
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SUBROUTINE  OUTPUT 

REAL  I EXACT,  IPHI,  K,  LENGTH,  MADMIT,  MEXACT,  MPHI , MVELP,  HYP, 
NORMAL 

COMPLEX  ADMIT,  ALPHA,  CEQN,  CEXACT,  CPHl,  CVEL,  CY,  IK.  1KSQ, 
TWOPIA,  VELP,  C INDATA 

COMMON  /I/  M,  NP,  NPP1 
COMMON  /R/  K,  PI,  TWOPI 
COMMON  /C/  ALPHA,  IK.  IKSO,  TWOPIA 
COMMON  /ID/  I CHECK  (102) 

COMMON  /RD/  LENCTII  (102),  NORMAL  (102,  2),  P (102,  2), 

. O < 102,  2.  2) 

COMMON  /CD/  CEQN  (102,  103),  CEXACT  (102),  CPHI  (102),  CVEL  (102), 

CY  ( 102) , CINDATA  ( 102) 

COMMON  /NOAUSS/  NCAUST 


DIMENSION  ADMIT  (102),  VELP  (102),  IEXACT  (102),  IPHI  (102), 

MADMIT  ( 102) , MEXACT  ( 102) . MPHI  ( 102) , MVELP  ( 102) , 
MYP  ( 102) , PADMIT  ( 102) , PEXACT  ( 102) , PPHI  ( 102) , 
PPVELP  ( 102) , PVELM  ( 102) , PVELP  ( 102) , PYP  ( 102) . 
REXACT  (102),  RPHI  (102) 


EQUIVALENCE 


(CEQN,  ADMIT) 
(CEQN(1,  3), 

( CEQN(  1 , 
(CEQN( 1 . 

( CEQN(  1 , 

( CEQN( 1 , 

( CEQN( 1, 

( CKQN( 1 , 

( CEQN ( 1, 


S)  , 
7)  , 
0)  , 
11) 
13)  , 
IS) 
17) 


IEXACT( 1)>, 
MADMIT(  1) ) . 
MPni(  D)  , 
MYP( 1>) , 
PEX\CT(  1 ) ) , 
PPVELP( 1) ) , 
PVELP( 1) ) , 
REXACTC 1 ) ) , 


( CEQN( I , 

2)  . 

<CEQN( 1. 

4)  , 

( CEQN< 1 , 

6)  , 

(CEQN( 1, 

B)  . 

(CEQN( 1, 

10) 

( CEQN( 1 . 

12) 

( CEQN( 1, 

14) 

(CEQN! 1, 

16) 

( CEQN( 1 , 

16) 

VELP(  D)  , 
IPHI(  D)  , 
MEXACT(  1 ) ) , 
MVELP ( 1)> , 
PADMIT( 1)) 
PPHI ( D)  , 
PVELM(  1 ) ) , 
PYP(  1) ) , 
RPHK  I)) 


PRINT  100 


FORMAT  ( "1",  6 IX, 

“**:*******  “ 

/ " ", 

6 IX. 

"* 

*"  / 

* ",  6 IX, 

”*  INPUT  *" 

/ " ", 

6 1 X, 

"* 

*"  / 

" ",  6 1 X, 

"********* 0 

////  " 

N" , 

15X, 

"EFFECTIVE  ADMITTANCE",  20X,  “ACOUSTIC  VELOCITY", 
"ACOUSTIC  POTENTIAL"  //  " “) 


22X, 


DO  1 I = 1,  NP 

REXACT(I)  = REAL  (CF.XACT(D) 

I EXACT!  I)  = A1MAG  (CEXACT!  I)) 

MEXACT(I)  = CABS  (CEXACT(D) 

PF.XACT(  I)  = 0.0 

IF  ( MEXACT!  I ) .NE.  0.0)  PEXACT(  I)  = ATAN2  (IEXACTO),  HEXACT(I)) 
MVELP ( I ) = CABS  ( CVEL(  I ) ) 

PVELP ( I)  =0.0 

IF  ( MVELP  ( I ) . NE.  0.0)  PVELP(I)  = ATAN2  (AIMAG  (CVEL(I)), 

REAL  ( CVEL(  I ) ) ) 

MYP!  1)  = CABS  (CY(  I) ) 

PYP( I)  = 0.0 

IF  (MYP!  I)  .NE.  0.0)  PYP(I)  = ATAN2  (AIMAG  (CY(D),  REAL  (CY(I))> 

IF  ( ICHECK(  I))  3,  4,  5 

CONTINUE 


WRITE  (6,  101)  I,  CEXACT(I) 

FORMAT  ( " ",  13,  B3X,  "(",  F13.  10,  V,  F13.10,  ■)  ■) 


CO  TO  22 
CONTINUE 


WRITE  (6,  102)  I,  CVEL(I) 

FORMAT  ( " ",  13,  40 X,  "(",  F13. 10.  *, *,  FI3.I0,  "> *> 


GO  TO  22 
CONTINUE 


WRITE  (6,  103)  I,  CY(  I) 

FORMAT  ( * ",  13,  9X,  *(",  F13. 10,  FI3. 10,  ’) ") 

CONTINUE 

CONTINUE 

NCAUST.  = 2 

NCAUST  = 2 * NCAUST 


WRITE  (6,  104)  K,  M,  ALPHA,  NGAUSZ,  NGAUST 


104  FORMAT  ( " I INPUT  FOR  THIS  CASE  IS: “ /// 

. ■ ",  25X,  “K=“,  F10.6,  5X,  “M=“,  15,  6X,  "ALPHA  = ( *, 

F5.  1,  F10.6,  ") * /-/ 

" ",  14X,  "NUMBER  OF  INTEGRATION  POINTS  IN  THE  XI  DIRECT 
“ ION  = ",  15,  // 

« “,  14X,  “NUMBER  OF  CAUSSIAN  POINTS  IN  THE  THETA  DIRECT 
“ION  = “,  15,  * (GAUSS  - LEGENDRE)”) 


PRINT  105 

105  FORMAT  ("1TIIE  CALCULATED  SURFACE  DISTRIBUTIONS  OF  THE  ACOUSTIC  P 

"OTENT1AL,  THE  ACOUSTIC  VELOCITY,  AND  THE  EFFECTIVE  ADM I 
"TTANCE  ARE: » //// 

'•  N",  11X,  “ P ( RHO , Z)“,  14X,  “PHI/COS(M*THETA) ",  20X, 

“VEL/COS(  M*THETA> “,  1BX,  “EFFECTIVE  ADM1TTNACE”  //  " ") 

DO  6 I = 1,  NP 
IF  ( I CHECK*  I) ) 7,  8,  9 
7 CONTINUE 

VELP( I)  = CPHI ( I) 
rVF.LM(I)  = CABS  (VELP(I)) 

PPVELP( I)  = O.0 

IF  (PVELM(I)  .NE.  0.0)  PPVELPC I)  = ATAN2  (AIMAG  (VELP(I)), 

REAL  ( VELP(  I)  ) ) 

CPHI(I)  = CEXACT(  I) 

RPni(I)  = REAL  (CPHKI)) 
irnui)  = aimag  (cpnid)) 

MPHKI)  = CABS  (CPHKI)) 

PPHK  I)  = 0.0 

IF  (MPHI(I)  .NE.  0.0)  PPHI(  I)  = ATAN2  (IPHKI),  RPHI(  I) ) 

ADMIT(I)  = VELP(  I)  /-CPHKI) 

IF  (MPni(I)  .EO.  0.0)  ADMIT(I)  = (0.0,  0.0) 

MADMIT(  I)  = CABS  (ADMIT(I)) 

PADMIT( I)  = 0.0 

IF  (MADMIT(I)  .NE.  0.0)  PADMIT(  I ) = ATAN2  (AIMAG  ( ADMIT(  D ) , 

REAL  ( ADMIT(  I) ) ) 

WRITE  (6,  106)  I,  P(I.  1).  P(I,  2),  VELP(I),  ADMIT(I) 

106  FORMAT  (“  ”,  13,  7X,  “(",  F5.3,  F6.3.  ")*,  44X,  2C13.5.  10X 

2C13.5) 


GO  TO  10 
8 CONTI NUE 

RPHKI)  = REAL  (CPHKI)) 

IPIIK  I)  = AIMAG  (CPHK  I)  ) 

MPHKI)  = CAES  (CPHKI)) 

PPHI(  I)  = 0.0 

IF  (NPni(I)  .NE.  0.0)  PPHK  I)  » ATAN2  (IPHKI),  RPHI<  I) ) 

VELP(I)  = CVEL( I) 

PVF.LM(  I ) = CABS  (VELP(I)) 

PPVELP( I)  = 0.0 

IF  (PVELMM)  .NE.  0.0)  PPVF.LP(  I)  = ATAN2  (AIMAG  (VELP(I)), 

REAL  ( VELP( l>) ) 

ADMIT(I)  = VELP(I)  / CPHKI) 

MADMIT(  I)  = GABS  (ADMIT(l)) 

PADMIT(I)  = 0.0 

IF  (MADMIT(I)  .NE.  0.0)  PADMIT(I)  = ATAN2  (AIMAG  (ADMIT(I))i 

REAL  ( ADMIT(  I))) 

WRITE  (6,  107)  I,  PM,  1),  P(  1 , 2),  CPHKI),  ADMIT(I) 

107  FORMAT  (*  ”,  13,  7X.  *(*,  F5.3,  F6.3,  ")",  8X.  2G13.5,  46X. 

2G13.5) 
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I 


GO  TO  10 
CONTINUE 

RPHKI)  = HEAL  < CPHK  I ) ) 

1PHI<  I)  = A I MAC  ( CPIIK  1)) 

Mrnu  n = cabs  <cphu  in 
min  i)  = o.o 

IF  (MTIIKI)  .NE.  0.0)  PPHI(I)  = AT  AN 2 (IPHI(I),  RPHI(I)) 

VELP1  I)  = CPHI ( I)  * CYC  I) 

PVF.LM1  I ) = CABS  (VELP(I)) 

PPVELPC I ) = 0.0 

IF  (PVELM1  I)  .NE.  O.O)  PPVELPC I)  = ATAN2  (AIMAG  CVELPCI)), 

REAL  (VELP( I))) 

ADM  I T(  I ) = CY(  I ) 

MADMIT(  l)  = CABS  (ADMIT(I)) 

P ADMIT! I)  = 0.0 

IF  ( MADMITC I ) .NE.  0.0)  PADMIT(I)  = ATAN2  (AIMAG  (ADMIT(I)), 

REAL  (ADMIT(  I))) 

WRITE  (6,  100)  I.  PCI,  1),  P(I,  2).  CPHKI),  VELP(I) 

FORMAT  ( “ ”,  13,  7X,  "(",  F5.3,  ",  ",  F6.3,  “)“,  8X,  2G13.5,  10X, 
2G13.5) 

CONTINUE 

CONTINUE 

PRINT  I 13 

FORMAT  ("ITHF.  MODULUS  OF  THE  ACOUSTIC  POTENTIAL,  THE  ACOUSTIC  VE", 
“LOCITY.  AND  THE  EFFECTIVE  ADMITTANCE  ARE:"  //// 

“ N”,  1 IX,  ”P(  RHO,  Z)“,  I4X,  ”PHI/COS(M*THETA) ",  20X, 

”VEL/COS(  M*THETA) " , I8X,  "EFFECTIVE  ADMITTANCE"  ✓ 

" ",  36X,  "EXACT  CALC",  20X,  "EXACT  CALC", 

20X.  "EXACT  CALC" 


DO  Id  I = 1.  NP 


IF  <ICIIECK<I>) 

17.  in. 

19 

17 

CONTINUE 

WRITE  (6, 

1 14) 

I,  P(  I, 

1)  , PC  I. 

- 

MYP<  I)  , 

MADMITt I) 

1 14 

FORMAT  < " 

",  13 

. 7X,  "( 

",  F5.3, 

. 2C13.5, 

10X,  2G13.5) 

GO  TO  20 

IB 

CONTINUE 

WRITE  (6, 

1 15) 

I.  P<  I, 

1)  , P<  I, 

■ 

MYP( I) , 

MADMIT(  I) 

1 15 

FORMAT  ( " 

",  13 

. 7X,  "( 

",  F5.3, 

G13.5, 

23X.  2G13.5) 

CO  TO  20 

19 

CONTINUE 

WRITE  (6. 

1 16) 

I,  M It 

1),  PC  I, 

• 

PVELM1  I ) 

. MYPC I) 

1 16 

FORMAT  ( " 

",  13 

1,  7X,  "( 

",  F5.3, 

2C13.5. 

10X.  GI3.5) 

20 

CONTINUE 

16 

CONTINUE 
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c 

c 

120 


C 

12 

C 


C 


C 


14 


0 


c 

15 
1 1 

C 


<; 


c 


c 

121 

C 


PRINT  120 


FORMAT  ( " 1THE  PRASE  OF  THE  ACOUSTIC  POTENTIAL,  THE  ACOUSTIC  VELO’, 
“CITY,  AND  THE  EFFECTIVE  ADMITTANCE  ARE:*  //// 

• N " , l IX,  ”P( RHO,  Z)“,  2 IX,  “PHI*.  33X,  "VEL",  24X. 
"EFFECTIVE  ADMITTANCE"  ✓ 

* ",  36X,  "EXACT  CALC",  20X,  "EXACT  CALC", 

20X,  "EXACT  CALC”  /✓  “ *) 


DO  11  I = 

1,  NP 

IF  ( I CHECK* I ) ) 

12,  13, 

14 

CONTINUE 

WRITE  (6, 

1 14) 

I,  P(  I. 

1) 

, P<  I 

PPVELP( 

n . 

PYPI 

GO  TO  15 
CONTINUE 

WRITE  (6, 

115) 

I.  PCI, 

i) 

. P<  I 

PYP(  I)  , 

P ADMIT* 

CO  TO  15 
CONTINUE 

WRITE  (6, 

1 16) 

I,  P<  I, 

1) 

. P(  I 

PPVELP(I),  PYP( 


, 2),  PEXACT(I),  PVELP(I), 
I),  PADMIT(I) 


, 2),  PEXACTC I)  . PPHI(I),  PVELP( I) , 

n 


, 2).  PEXACT(I),  PPHI(I),  PVELP(I), 
I) 


CONTINUE 

CONTINUE 


WRITE  THE  SURFACE  DISTRIBUTIONS  TO  THE  OUTPUT  FILE. 


WRITE  (10.  121)  (CPHI(I).  VELP(I).  I = 1.  NP> 

FORMAT  (4020. 10) 


RETURN 

END 
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The  following  computer  program  reads  the  previous  programs 
output  file  containing  the  values  of  the  acoustic  quantities  on  the  surface  of 
the  body  and  calculates  the  values  of  the  acoustic  quantities  at  any  given 
point  in  the  field  surrounding  the  body.  The  required  inputs  are  the  points 
where  the  acoustic  quantities  are  known  on  the  surface  of  the  body  Q,  the 
normals  to  the  body  at  these  points  NQ,  the  length  of  each  integration 
interval  LENGTH,  the  points  in  the  field  where  the  acoustic  quantities  are 
required  P,  and  some  arbitrary  normal  at  these  points  NP  as  the  normal 
acoustic  velocity  is  calculated.  The  problem  specification  data  is  also 
required  again;  that  is,  k the  wave  number  and  m the  mode  number. 

The  program  prints  out  all  the  geometric  input  data  and  the 
acoustic  potential  and  normal  acoustic  velocity  at  the  field  points.  It  also 
calculates  and  prints  out  the  SPL(dB)  at  each  field  point. 


1 


PHQCRAM  EXPFF  (INPUT.  OUTPUT.  TAPE 10,  TAFEll, 
TAPES  = INPUT.  TAPE6  = OUTPUT) 


G 

C 


c*********************************************************************** 


c* 

c* 

c* 

c* 

c* 

c* 

c* 


THIS  PROGRAM  CALCULATES  THE  ACOUSTIC  POTENTIAL  AND  THE 
ACOUSTIC  VELOCITY  IN  THE  FIELD  SURROUNDINC  ANY  AX I SYMMETRIC 
BODY  EMPLOYING  THE  SURFACE  DISTRIBUTIONS  OF  THE  ACOUSTIC 
POTENTIAL  AND  TIIF.  NORMAL  ACOUSTIC  VELOCITY. 


* 

* 

* 

* 

* 

* 

« 


C************** ****** *************************************************** 
C*  * 

C*  A CYLINDRICAL  FORMULATION  OF  THE  PROBLEM  IS  EMPLOYED.  * 

C*  * 

c*  * 

c****************»****************************************************** 

c 

c 

REAL  K.  LENGTH 

COMPLEX  CPU  IP,  CPHIQ.  CVELP.  CVELQ.  IK.  IK5Q 


COMMON  /I/  M.  NP,  NQ 

COMMON  /R/  FOURPI.  K,  PI.  TVOP I 

COMMON  /C/  IK.  1KSQ 

COMMON  /RD/  LENCTH  (102).  P (0,  2).  PNORMAL  (9,  2), 

Q (102.  2).  QNORMA1.  (102,  2) 

COMMON  ✓CD/'  CPHIP  (9).  CPHIO  (102),  CVELP  (9).  CVELO  (102) 
COMMON  /NGAUSS/  NGAUSS 
COMMON  /GAUSS/  GAUSST  (48,  2) 


C 

c 

c 


CALL  INPUT 
REAL  CONSTANTS 


TWOPI  = 20  * PI 
FOURPI  - 4.0  * PI 


C 

C 

C 


COMPLEX  CONSTANTS 

* K 


IK  = (0.0.  1.0) 
IKSU  = IK  * IK 


CALL  CALC 
CALL  OUTPUT 
STOP  "NORMAL" 


END 
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BLOCK  DATA 


REAL  K 


COMMON  /\/  M,  NP.  NO 
COMMON  /R/  FOURPI.  K,  PI,  TWOPI 
COMMON  /NCAUSS  ' NGAUSS 
COMMON  /GAUSS/  CAUSST  (48,  2) 

THE  NUMBER  OF  POINTS  IN  THE  FIELD  WHERE  THE  ACOUSTIC  POTENTIAL 
AND  THE  ACOUSTIC  VELOCITY  ARE  TO  BE  CALCULATED. 


DATA  NP  / 0 / 

THE  NUMBER  OF  POINTS  ON  THE  SURFACE  OF  THE  BODY  WHERE  THE 
ACOUSTIC  POTENTIAL  AND  THE  NORMAL  ACOUSTIC  VELOCITY  ARE  KNOWN. 


DATA  NQ  / 102  / 

DATA  PI  / 3.1415926535898  X 
DATA  NCAUSS  / 48  / 


DATA 


DATA 


DATA 


DATA 


( (CAUSST( I , J)  , J = 
/ 0.01627674484960, 
0.04881298513605, 
0.081 29749546443 . 
0. 1 136958501 1067, 
0.  14597371465490, 
O. 17809688236762. 
0.21003131046057, 
0.24174315616384. 
0.27319031259105. 
O . 30436494433450 , 
O . 33520852289263 , 
0.36569686147231. 
((CAUSSTt  I . J)  . J = 
/ 0.39579764982891, 
0 . 42547898840730 , 
0.45470942216774, 
0 . 48345797392060 , 
0.5  I 169417715467, 
0. 539303 10832436 , 
0.56651041856140, 
0 . 59303236477757 , 
0.61892504012547, 
0.6441 6340378497 , 
0.66071831004392. 
0.69256453664217, 
< (CAUSSTC I . J) . J = 
/ 0.71567681234897, 
O . 73803064374440 , 
0.759602341 17665, 
0.78036904386743, 
0.80030874413914, 
0.81940031073793, 
0.83762351 122819, 
0 . 85495903343460 , 
0.87I38350590930, 
0.88609451740242, 
0.90  1 4606353  1 585. 
0.91507142312090, 
((GAUSSTO,  J),  J = 
/ 0.92771245672231, 
0.93937033975276, 
0.95003271778444, 
O. 95968829 1 44874, 
O . 06332682046326 . 
0.97593917458514, 
0.08251726356301 , 
0.908054 12632962, 
O . 99254390032376 , 
0.99590104298721 , 
0.09036437586310, 
0 . 09968930308323 , 


1,  2),  I = 1,  12) 

0 . 0325506 1449236 , 
0.032516I1871387, 
0.03244716371406, 
0.03234382256858, 

0 . 03220620479403 , 
0.03203445623199, 
0.03182875889441 , 
0.03158933077073, 
0.03131642559686, 
0.03101033258631 . 
0.03067137612367, 

0. 0302999  1542083  / 

1.  2).  I = 13.  24) 
0.02989634413633, 
0.02946108995817, 
0.02399461415056, 

0 . 02H4974 1 106509, 
0. O27970007616B5, 
0.02741296272603, 
0.02682686672559, 
0.02621234073567, 
0.02557003600535, 

0 . 02490063322248 , 
0.02420484179236, 

0. 02348339908593  / 

1,  2) , I = 25,  36) 
0.02273706965833. 
0.02196664443874, 
0.02117293989219, 
0.02035679715433, 
0.01951908114013, 
0.01866067962741, 
0.O17782502316O5, 
0.01688547986425, 
0.01597056290236, 
0.01303872102699, 
0.01409094177231 , 

0. 01312822956696  / 

1,  2).  I = 37,  48) 
0.01215160467109. 
0.011 16210209984. 
0.01016077053501 , 
0.00914867123078, 

O . 003 1 2687692570. 
0.00709647079115, 
0.00603854550424, 
0.00301420274293. 

0 . 00396455433844 , 
0.00291073181793, 
0.00185396078893, 
0.00079679206555  / 


END 
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SUBROUTINE  INPUT 
REAL  K,  LENGTH 

COMPLEX  CPHIP , CPHIO,  CVEI.P.  CVELQ 

COMMON  /I/  M,  NP.  NO 

COMMON  /R/  EOURPI.  K,  PI,  TWOPI 

COMMON  /RD/  LENGTH  (102),  P (9,  2),  PNORMAL  (9,  2), 

U (102,  2),  QNORMAL  (102,  2) 

COMMON  /CD/  CPHIP  (9),  CPHIQ  (102),  CVELP  (9),  CVELQ  (102) 

READ  (5,  100)  ((P(I,  J) , J = 1,  2),  I = 1,  NP) 

READ  (5.  100)  ((PNORMAL  (I,  J) , J = 1 , 2)  , I = 1.  NP) 

READ  (5.  100)  ((Q( I.  J) , J = 1,  2) , I = 1 , NQ) 

READ  (5.  100)  ((QNORMAL  (I,  J) , J = 1,  2),  I = 1,  NQ) 

READ  (5,  100)  (LENGTH(I),  1=1,  NQ) 

FORMAT  ( BG 10.0) 

PRINT  101 

FORMAT  ( " 1GE0MF.TR I C INPUT  DATA:  " /// 

» ",  4X,  "N\  2 IX,  "Q(RHO.  Z>",  38X,  "NORMAL!  RBO,  Z)  * , 

29X.  “LENGTH"  //  " ") 

WRITE  (6,  102)  (I.  Q(  1 , 1),  Q(  1 , 2),  QNORMAL! I , I),  QNORMAL! I,  2), 
LENGTH!  I),  I = 1 , NQ) 

FORMAT  ("  ",  2X,  13,  10X,  "<".  F13.  10,  F13.10,  ')",  21X. 

•(",  F13.10,  F 13. 10 , 17X, 

. F13.10) 

PRINT  103 

FORMAT  ( “ 1 FIELD  POINT  INPUT  DATA:  " /// 

“ ",  4X,  "N",  21X,  "P(RHO,  Z) ",  38X,  "NORMAL!  RHO,  Z) " // 

n » j 

WRITE  (6,  104)  (I,  P(I,  1),  P(I,  2).  PNORMAL!  1 , 1),  PNORMAL!  I,  2), 
I = 1,  NP) 

FORMAT  ( " ",  2X.  13,  9X,  "(",  F14.10,  F14.10,  ")",  20X. 

"(  ",  F 13 . 10.  ",  ",  F 13 . 10.  ")  ") 

READ  ( lO,  105)  (CPHIQ(I),  CVELQ(I),  I = 1 , NQ) 

FORMAT  (4G20.10) 

PRINT  106 

FORMAT  ("t“,  6 IX,  "********"  / * *,  61X,  "*  **  / 

. " ",  6 1 X,  "*  BODY  61X,  "*  *•  / 

» ",  6 IX,  "********"  ////  " N",  2 IX,  ”Q( RHO,  Z) * , 25X, 

“ACOUSTIC  POTENTIAL",  2 IX,  "ACOUSTIC  VELOCITY"  //  " ") 

WRITE  (6.  107)  (I,  Q(l,  1),  Q(  1 , 2),  CPHIQ(I),  CVELQ(  I) , 

I = 1,  NQ) 

FORMAT  ( “ ",  13.  9X,  "(",  F13.10,  F13.10,  *)", 

9X.  *(",  FI4.10.  F14.10,  *)“, 

8X.  "(",  F14.10,  F 14. 10 , *) “) 


READ  (II,  1 100)  M,  K 

FORMAT  (120,  G20.0) 

RETURN 

END 
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SUBROUTINE  CALC 
C 

REAL  K,  LENCTH,  NDOTN,  NRHOQP.  NZPZD,  RZQZD,  NZQP 

COMPLEX  CPniP,  CPHIQ,  CVELP,  CVELO,  C,  GP,  GPP,  IK,  IKSQ,  II,  12, 
. 13,  14 

C 

COMMON  /I/  M,  NP.  NQ 

COMMON  /R/  FOURPI,  K,  PI.  TWOPI 

COMMON  /C/  IK,  IKSO 

COMMON  /RD/  LENGTH  (102),  P (9,  2).  PNORMAL  (9,  2), 

Q (102,  2),  QNORMAL  (102,  2) 

COMMON  /CD/  CPHIP  (9).  CPHIQ  (102),  CVELP  (9),  CVELQ  (102) 

COMMON  /NCAUSS/  NCAUSS 
COMMON  /GAUSS/  GAUSST  (48,  2) 

C 

DO  1 I = 1.  NP 
CPHIP(I)  = (0.0.  0.0) 

CVELP ( I)  = (0.6,  0.0) 

XI  INTEGRATION 

DO  2 J = I,  NQ 
GAUSZ  = TWOPI  * T.F  GTH(J) 

ZD  = Q(J,  2)  - P(  2) 

NZQZD  = QNORMAL(J,  2)  * ZD 
NZPZD  = P NORMAL!  1 , 2)  * ZD 
ZSQ  = ZD  * ZD 

IU10SQ  = (Q(J.  1)  - P(I.  1 ) ) **2 
RH0QP2  = 2.0  * Q(J,  l)  * P(I.  1) 

NRHOQP  = QNORMAL(  J , I ) * PN0RMAL(  I , I ) 

NZQP  = QNORMAL(  J , 2)  * PNORMAL!  1 , 2) 

UHOZSQ  = RHOSQ  + ZSQ 

THETA  INTEGRATION 

DO  3 IT  = 1 . NCAUSS 
THETA  = IM  * GAUSST! IT.  1) 

GAUSZT  = GAUSZ  * GAUSST!  IT.  2) 

C 

COST  = COS  (THETA) 

COSMT  = COS  (M  * THETA) 

C 

R = SORT  ( RUOZSQ  ♦ R1IOQP2  * (1.0  - COST) ) 

DRDNP  = ( PNORMAL(  1 , 1)  * (P(I,  I)  - Q(J,  1)  * COST)  - NZPZD)  / R 
DRDNQ  = ( ONORM\L(  J . 1)  * ( Q( J . 1)  - P( 1 , 1)  * COST)  4 NZQZD)  / R 
NDOTN  = NRHOQP  * COST  + NZQP 
C 

G = Q(J.  I)  * GAUSZT  * COSMT  * CEXP  ( IK  * R)  / R 
OP  ■ G * (IK-  (1.0  / R)) 

CPP  = G * ( 1KSQ  - (3.0  * IK/  R)  + (3.0  / (R  * R))) 

C 

11  = GP  * DRDNQ 

12  * -C 


13  = GPP 

* 

DRDNQ  * 

DRDNP  - 

(GP  * NDOTN 

/ 

R) 

14  = -CP 

* 

DRDNP 

CPHIP(  I) 

S 

CPHIP!  I) 

I 1 * 

CPHIQ!  J)  + 

12 

* 

CVELQ(J) 

CVF.LPt  I) 

s 

CVELP(  I) 

■f 

13  * 

CPHIQ(J)  + 

14 

* 

CVELQ(  J) 

CONTINUE 

CONTINUE 

CPHIP( I) 

U 

cpniP(  i) 

/ 

FOURPI 

CVELP!  I ) 

= 

CVELP  ( I) 

/ 

"OURPI 

C 

1 CONTINUE 

C 

RETURN 

END 
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■ 


c 

c 


c 

c 


c 

c 


SUBROUTINE  OUTPUT 


REAL  I PHI , IPHIP.  IVEL,  IVELP,  K,  LENGTH,  MPHI.  HPHIP,  HVEL, 
MVELP 

COMPLEX  CPHIP,  CPHIQ,  CVELP,  CVELQ,  Y 


COMMON  /I/  M,  NP.  NQ 

COMMON  /R/  FOURPI,  K,  PI,  TWOPI 

COMMON  /RD/  LENGTH  (102),  P (9,  2),  PNORMAL  (9,  2), 

Q (102,  2),  QNORMAL  (102,  2) 

COMMON  /CIV  CPHIP  (9).  CPHIQ  (102),  CVELP  (9),  CVELQ  (102) 
COMMON  /NGAUSS/  NGAUSS 


DIMENSION  I PH  I (9),  IPHIP  (9).  MPHI  (9),  MPHIP  (9). 

PPHI  (9),  PPHIP  (9).  RPHI  (9),  RPHIP  (9), 
IVEL  (9),  IVELP  (9),  MVEL  (9),  MVELP  (9), 
PVEL  (9),  PVELP  (9),  RVEL  (9),  RVELP  (9) 


INITIALIZE  EXACT  SOLUTION. 


1 , NP 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


DO  2 I 
RPUI ( I)  = 

IPHI(I)  = 

MPHK  I)  = 

PPIIK  I)  = 

RVEL(  I)  = 

I VF.L(  I ) = 

MVEL(  I)  = 

PVF.L(  I ) = 

CONTINUE 

NGAUSS  - *>.  * NGAUSS 


WRITE  (6,  100)  K,  M,  NCAUSS 


100 


/// 


FORMAT  ( "1  INPUT  FOR  THIS  CASE  IS: 

» ",  60X,  "K  = “,  F10.6  //  " ”,  60X,  "M  =",  15  // 

“ “,  1 3X,  "NUMBER  OF  GAUSSIAN  POINTS  IN  THE  THETA  DIRECT" 

“ION  = ",  15,  » (GAUSS  - LEGENDRE)") 


= 1 


DO  5 
nPHIP( I) 
RVELP  ( I ) 
IPniP( I) 

I VELP( I) 

MPH I P(  I) 
MVEI.P(  I ) 
PPHIP(  I) 
PVELP  ( I ) 

IF  ( MPIIIPI  I) 
IF  ( MVELP*  I) 
CONTINUE 


NP 

REAL  ( CPHIP( I)) 

REAL  ( CVELP ( I) ) 

A I MAG  ( CPHIP(  D) 

A I MAG  ( CVELP ( I)) 

CABS  (CPIIIP(D) 

CABS  ( CVELP ( I ) ) 

0.0 

0.0 

.NE.  0.0)  PPHIP( I) 
. NE.  0.0)  PVELP(I) 


= AT AN 2 
= AT AN 2 


( IPHIP( I) , 
( IVELP( I) , 


RPHIP( I)) 
RVELP ( I) ) 


PRINT  103 


103  FORMAT  ( ”1  ' 


53X,  "**************"  / ’ ", 
58X,  "*  CALCULATED 

5BX,  “**************"  ////  ' 
25X,  "ACOUSTIC  POTENTIAL", 

2 1 X,  "ACOUSTIC  VELOCITY"  // 


38X,  "* 
58X.  "* 
N", 


*"  / 
*•  / 

21X,  *P( RHO,  Z)", 


■) 


102 


WRITE  (6. 

102)  (I,  P( 

I. 

1),  P( 

I. 

2).  CPHIP( 

I) 

• 

I = 1 

* 

NP) 

FORMAT  ( * 

",  13,  8X, 

"( 

",  F14. 

10, 

, F14. 

10 

, 

8X, 

"( 

",  F 14 . 

10. 

, F 14 . 

10 

, 

BX, 

*( 

",  F14. 

10, 

F14. 

II 

•>  *, 

")  ", 


MM 


J 
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r.i.m  109 

FORMAT  ( " 1THE  MODULUS  OF  THE  ACOUSTIC  POTENTIAL  AND  THE  ACOUSTIC* , 
" VELOCITY  ARE:"  ////  * N “ , 1 IX,  "PCRHO,  Z)  * , 14X, 
"Pni/COSC M*THETA) ",  20X,  "VEL/COSC  M*THETA) “ , 8X, 

"SPL  ( PB) " , 16X.  “Y"  / 

• ",  36X,  "EXACT  CALC",  20X,  "EXACT  CALC"  // 

" * ) 

DO  7 I = l,  NP 

SPL  = 20.0  * ALOG10  (K  * MPHIPCI))  + 146.6 
Y = CVKLP( I)  / CPHIP( I) 

WRITE  (6,  110)  I,  P(I,  1),  PH.  2),  MPHICI),  MPHIPCI),  MVEL(I), 
MVELP(I),  SPL,  Y 

FORMAT  ("  ",  13,  6X,  F6.3,  F7.3,  7X,  2G13.5,  10X, 

2G13.5,  2X.  F7.2,  5X,  2G13.5) 

CORTINUE 
PRINT  112 

FORMAT  ( "1THE  PHASE  OF  THE  ACOUSTIC  POTENTIAL  AND  THE  ACOUSTIC  V*. 
"ELOCITY  ARE:"  ////  * N",  11X,  ”P(RHO,  Z) “ , 21X, 

"rni",  33X.  "VEL"  / 

“ ",  36X,  "EXACT  CALC",  20X,  "EXACT  CALC"  // 

h » y 

WRITE  (6,  113)  (I,  PCI,  1),  PCI,  2),  PPHICI),  PPHIPCI),  PVELC I) , 
PVELPCI),  I 1 1,  NP) 

FORMAT  C"  ",  13.  6X.  "C",  F6.3,  *,*,  F7.3,  ")*,  7X,  2C13.5,  10X, 
2GI3.5) 

RETURN 

END 


i 
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Paper  accepted  for  publication  in  the  Journal  of  the  Acoustical 
Society  of  America 

"Prediction  of  the  Sound  Field  Radiated  From  Axisymmetric 
Surfaces" 


PREDICTION  OF  THE  SOUND  FIELD  RADIATED  FROM 
AXISYMMETRIC  SURFACES 


W.  L . Meyer*  W.  A.  Bell',  M.  P.  Stal lybrass^and  B.  T.  Zinin" 

*School  of  Aerospace  Engineering 
**School  of  Mathematics 

Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 

Abs  tract 

A general  analytical  method  for  determining  the  radiated  sound  fields 
from  axisymmetric  surfaces  of  arbitrary  cross  section  with  general  boundary 
conditions  is  developed.  The  method  is  based  on  an  integral  representation 
for  the  externa  1 solutions  of  the  Helmholtz  equation.  An  integral  equation 
is  developed  governing  the  surface  potential  distribution  which  gives  unique 
solutions  at  all  wave  numbers.  The  axisymmetric  formulation  of  the  problem 
reduces  its  solution  to  the  numerical  evaluation  of  line  integrals  by 
Gaussian  quadrature.  The  applicability  of  the  solution  approach  for  both  a 
sphere  and  finite  cylinder  is  demonstrated  by  comparing  the  numerical  results 
with  exact  analytical  solutions  for  both  discontinuous  and  continuous  boundary 
conditions.  The  method  is  then  applied  to  a jet  engine  inlet  configuration 
and  the  computed  results  are  in  good  agreement  with  exact  values. 

I . Introduction 

To  reduce  the  noise  radiated  to  the  community  from  turbofan  inlets,  the 
effects  of  sound  suppression  material  in  the  inlet  and  the  spatial  distri- 
bution of  the  sound  source  on  the  radiated  sound  levels  and  patterns  must  be 
determined.  Analytical  techniques  for  predicting  these  effects  must  be  cap- 
able of  dealing  with  general  axisymmetric  geometries  and  complicated  boundary 
conditions  which  are  encountered  in  multiply- lined  inlets.  For  instance,  in 
a typical  inlet  the  compressor- fan  combination  represents  a noise  source 
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with  a nonuniform  spatial  excitation  pattern.  Thus,  the  analytical  method 
should  be  capable  of  taking  into  account  sound  sources  of  general  spatial 
distribution.  Also,  inlets  may  contain  multiple  acoustic  liners  to  reduce 
the  radiated  sound  power  and  admittance  boundary  conditions  are  commonly 
used  to  account  for  the  absorption  characteristics  of  the  liner.  Therefore, 
the  analytical  method  must  be  capable  of  dealing  with  spatially  varying  sur- 
face admittances.  Finally,  the  method  should  be  capable  of  predicting  the 
characteristics  of  the  radiated  sound  field  in  an  infinite  domain.  Keeping 
these  requirements  in  mind,  the  work  presented  in  this  paper  describes  the 
results  of  an  investigation  which  has  been  concerned  with  the  analytical 
determination  of  radiated  sound  fields  from  axisymmetri:  surfaces  of  arbi- 
trary cross  section  and  with  general  boundary  conditions. 

The  method  used  in  this  investigation  is  based  on  an  integral  form  of 

1-6 

the  solutions  of  the  Helmholtz  equation.  With  this  formulation  the  acous- 
tic potential  anywhere  external  to  the  surface  can  be  found  once  the  distri- 
bution on  the  surface  is  known.  Thus,  to  determine  the  radiated  sound  field 
the  problem  reduces  to  the  determination  of  the  distribution  of  the  acoustic 
potential  on  the  two  dimensional  surface  of  the  geometry  under  consideration 
instead  of  solving  the  Helmholtz  equation  in  the  surrounding  infinite  three 
dimensional  domain. 

It  has  been  previously  shown  that  when  applied  to  exterior  sound 
radiation  problems  the  classical  techniques  fail  to  produce  unique  solutions 
at  frequencies  corresponding  to  certain  interior  eigenvalues  of  the  geometries 
under  consideration.  Unless  special  precautions  are  taken,  straight-forward 
numerical  solution  of  the  integral  equation  produces  large  errors  at  frequen- 
cies close  to  these  eigenvalues.  For  the  general  geometries  of  interest  in 
this  study,  these  eigenfrequencies  are  not  known  a priori.  Therefore, 


I 
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the  frequencies  about  which  large  numerical  errors  can  occur  cannot  be  easily 

avoided.  A critical  review  of  available  analytical  techniques  for  avoiding 

these  errors  is  provided  by  Burton  in  Ref.  1.  In  a search  for  an  appropriate 

technique  for  use  in  the  present  study  of  inlets,  the  authors  programmed  each 

of  these  methods  for  a sphere  and  obtained  numerical  results  for  the  surface 

and  radiated  sound  field.  This  study  showed  that  the  method  of  Burton  and 
4 

Miller  was  the  most  straightforward  to  implement.  However,  an  inter- 
pretation of  a strongly  singular  integral,  given  in  the  analysis  in  Ref.  5 
by  Meyer,  et.al.  was  necessary  for  the  equations  to  be  amenable  to  numerical 
solution.  Basically  the  method  proposed  by  Burton  and  Miller  involves  a 
reformulation  of  the  integral  equation  for  the  acoustic  potential  and  the 
solutions  obtained  are  valid  at  all  frequencies.  It  also  yields  the  most 
consistently  accurate  results  for  a given  number  of  points  at  which  the 
acoustic  potential  is  numerically  evaluated  on  the  surface.  Therefore,  the 
method  based  on  the  analysis  in  Ref.  5 has  been  chosen  for  this  investigation. 

The  resulting  integral  equation  for  the  surface  acoustic  potential  is 
solved  numerically  and,  for  axisymmetric  geometries,  the  equation  reduces  to 
the  evaluation  of  line  integrals.  Thus,  the  axisymmetric  case  can  be  reduced 
to  an  equivalent  one-dimensional  problem.  Having  discretized  the  integral 
equation,  the  resulting  system  of  algebraic  equations  is  solved  using  complex 
Gauss- Jordan  elimination.  Since  the  coefficient  matrix  involves  the  free 
space  Green's  function,  which  becomes  singular  as  two  points  on  the  surface 

approach  one  another,  numerical  techniques  are  presented  which  can  deal 
with  these  singularities  and  yield  accurate  results.  Gaussian  integration  is 
used  to  increase  the  accuracy  of  the  solution  without  significant  penalties 
in  computer  storage  and  time  requirements.  The  applicability  of  the  integral 


formulation  and  the  accuracy  of  the  numerical  techniques  are  demonstrated  by 
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computing  the  surface  and  far  field  distributions  of  the  acoustic  potential 
on  both  a sphere  and  a finite  cylinder.  The  numerical  results  are  compared 
with  known  exact  solutions  generated  by  the  separation  of  variables  technique. 
Surfaces  with  spatially  varying  forcing  functions  and  admittances  are 
considered,  for  different  tangential  modes,  to  evaluate  the  capability  of  the 
integral  approach  to  handle  boundary  conditions  of  a general  nature.  With 
the  sphere,  agreement  between  computed  and  exact  results  is  to  three  signifi- 
cant figures.  For  the  cylinder  agreement  is  to  two  significant  figures.  The 
effect  on  the  accuracy  of  discontinuous  boundary  conditions  involving  nonzero 
admittances  over  the  surface  and  of  the.  corners  encountered  in  the  cylindrical 
configuration  are  also  presented.  Finally,  the  numerical  results  for  an 
inlet  configuration  are  compared  with  exact  solutions  and  agreement  is  to 
within  ten  per  cent. 


II . Theoretical  Considerations 

In  this  section  the  general  three  dimensional  integral  representation  of 
the  solutions  of  the  Helmholtz  equation  is  developed  for  application  to  radi- 
ation problems.  This  particular  formulation  yields  unique  solutions  at  all 

frequencies  and  does  not  have  strong  singularities  which  are  difficult  to 
handle  numerically.  The  general  integral  equation  is  then  specialized  for 
axisymmetric  geometries.  A more  detailed  development  is  presented  in  Ref.  5. 
General  Formulation 

Beginning  with  the  three  dimensional  Helmholtz  equation  which  governs 
the  spatial  dependence  of  the  acoustic  field  for  harmonic  oscillations 

2 2 m 
ycp+kcp=0  v ' 


X 


5 


where  cp  is  the  acoustic  potential  and  k is  the  wave  number;  the  standard 


integral  representation  of  the  exterior  potential  is  found  to  be 
J J (C?(Q)  " G(P»Q)  )dSq  = 4n'P(P) 


1,6 


(2) 


where  the  term  — represents  an  outward  normal  derivative  with  respect 

sn, 

to  the  body  S as  shown  in  Fig.  1;  that  is. 


^ *«>  • \ 
q 


(3) 


Also,  G(P,Q)  is  a fundamental  three  dimensional  solution  of  the  Helmholtz 
equation  and  is  taken  to  be  the  free  space  Green's  Function  for  a point 


source0  defined  as 


G(F,Q)  = 


,ikr(P,Q) 
r (P,Q) 


(4) 


From  Eq,  (2),  if  the  acoustic  potential  and  the  normal  acoustic  velocity 


are  known  at  each  point  on  the  surface  of  the  body  then  the  acoustic 


potential  may  be  calculated  anywhere  in  the  exterior  domain. 

To  solve  for  the  surface  potential,  the  point  P is  moved  to  the  surface 
of  the  body,  and  Eq.  (2)  then  becomes 

I J'  <»«>  ' G<p-Q)  ^ )dsq  - 2 


(5) 


For  the  inhomogeneous  Robin  boundary  condition  employed  in  this  study, 
a relation  between  dcp(Q)/dnq  and  cp (Q)  exists  which  is  given  by 


- *«)<?«)  = A(Q), 

q 


(6) 


so  that  Eq.  (5)  can  be  written  in  terms  of  the  potential  only;  that  is, 

J j *«>  ds,  - J j (7) 


2 tt  cp(P)  + j J A(Q)  G(P,Q)  dSc 


If  the  acoustic  velocity  A(Q)  and  the  admittance  Y ( Q ) are  specified  at  each 
point  on  the  surface  of  the  body,  then  the  acoustic  potential  may  be  calcu- 
lated at  each  point  using  Eq.(7). 

As  mentioned  earlier  this  equation  does  not  yield  unique  solutions  when 

the  wave  number  k is  an  internal  eigenvalue  associated  with  the  problem  under 

consideration.  Since  these  eigenvalues  are  not  known  a priori  for  general 

bodies,  the  formulation  cannot  be  relied  upon  to  give  consistently  good  results. 

2 3 4 

There  are  a number  of  papers  in  the  literature  ’ ’ dealing  with  this  problem, 
and  the  relative  merits  and  shortcomings  of  the  methods  employed  are  discussed 
in  detail  in  Ref . 1 . 

An  attractive  approach  from  an  analytical  point  of  view  is  provided  by 

4 ... 

Burton  and  Miller  who  have  suggested  the  use  of  the  following  identity  to 

derive  an  alternative  integral  equation  for  the  acoustic  potential  at  the 


surface . 


’SP  - J J[»< 


Bn  Bn 
p q 


MSI  -ld 

Bn  J 


This  equation  can  now  be  solved  for  cp(P)  by  using  Eq.  (6)  to  relate  the  normal 
acoustic  velocity  and  the  potential  at  the  surface.  However  this  integral 
equation  has  its  own  set  of  associated  eigenvalues  at  which  unique  solutions 
cannot  be  obtained.  To  circumvent  the  problem  associated  with  the  solution  of 
the  integral  equations  derived  from  Eqs.(5)and  (8),  Burton  and  Miller  suggested 
the  solution  of  the  following  linear  combination  of  these  equations: 


J J 


- G(P,Q) 


MSI 
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- « J J (cp(Q)  ) dS 

J J on  dn  -n  in  q 


p q 


2n  (c?(P)  + a ) 

on  p 


(9) 


whore  3 qp  / 9n  and  Ware  related  by  Eq . (5).  Equation  (9)  will  yield  unique 
solutions  if  the  complex  coupling  constant  is  properly  chosen.  It  is  shown 

that  pi  must  moot  the  following  restrictions  to  guarantee  that  Eq . (9)  yields 

. * 
unique  solutions: 


Im  (cc)  ^ 0 

k real  or  Imaginary 

Im  ( a ) = 0 

k complex 

(10) 


A problem  arises  in  the  numerical  solution  of  Eq.  (9)  as  the  third  term 

on  the  left  hand  side  is  strongly  singular  ir.  its  present  form  as  the  point 

Q approaches  the  point  P on  the  surface  of  the  body.  The  authors  of  this 

paper  have  shown  that  this  difficulty  can  be  overcome  by  a proper  interpre- 

5 8 

tation  of  this  singular  term.  Employing  a vector  transformation  and  taking 
the  Cauchy  Principle  Value,  Eq.  (9)  is  shown  to  be  equivalent  to 

J J [*«>  ^ ^ ] ", 

s q q 

+ O J J [*«>  - *W]  §^j;1<a  -JS 

- or  <p(P)  J J(np  . nq)(ik)2  G(P,Q)  dSq 
S 


It  has  been  pointed  out  to  us  by  a reviewer  that  an  equation  of  the  same 
general  form  as  Eq.  (9)  has  been  given  by  Chertock  . However,  for  an  arbi- 
trary, smooth  surface,  Chertock  did  not  interpret  this  integral  equation 
correctly.  Specifically,  the  limit  indicated  in  the  final  term  of  (A17), 
Ref.  (7),  does  not  exist  as  may  be  verified  for  the  simple  case  of  a sphere. 
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j 


"“I  J 


dS  = 2 n [cp(P)  + a 1 

dn  q LV  dn  J 


All  of  the  terms  in  Eq.  (11)  are  now  well  defined;  however,  all  the  inte- 
grands are  oscillatory  and  singular  so  that  care  must  be  taken  in  their 
numerical  approximation, 

Axisymmetric  Formulation 

When  dealing  with  a body  of  revolution  as  sho\n  in  Fig.  2 an  axisym- 

9 

metric  formulation  of  the  problem  is  advantageous.  This  being  the  case  an 
element  of  area  dSq  becomes  pdsd0  where  s is  the  distance  along  the  perimeter 
of  the  surface  in  the  p-z  plane.  Assuming  an  acoustic  velocity  distribution 


of  the  form 


v(s)  cos  m 0 


and  describing  the  s dependence  of  the  potential  function  by 


$ (s)  3 


cos  m Q 


and  letting  0=0  (so  that  cos  0 = 1)  Eq.  (11)  becomes: 


J J *<V 


cos  m 0 dS 

dnq  q q 


a *(sp)  j J G(P,Q)(ik)2(np.nq)  dSq 
S 

2 

+ « J J l4(V  cos » 0,  - #<sP>]  trj 


- J I »<v  G(PiQ)  cos  m Gq  dS^ 
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- a [ f v(s  ) — cos  m 0 dS 

J J v q'  dn  q q 

S p 

= 2r[_^(sp)  + O'  v(s^)J 


Now,  three  sets  of  functions  are  defined: 
Influence  Functions 


I,  (r  ) = 2 j G(P,Q)  cos  m 0 dQ 
x F4  v o 4 4 

I2  (r  ) = 2 a J cos  m 0 d0 

L pq  J on  q a 


o p 


Kernel  Functions 


K,  (r  ) = 2 [ 

1 p q i 


^P-».Q)  cos  tn  0 dQ 

b Snq  q q 


K?  (r  ) = 2 a [ 

2 pq  Jo 

Forcing  Functions 

rn  2 

F1  (rpq>  = 2a  J G<P’Q)<ik)  <VV 


cos  m 0 d0  0 t e 
bnpsnq  q q’  q p 


dG„ 


F (r  ) » 2«  [ d 0 , 0^0 

2 pq  bn  an  q °q  r ui 


(15) 


(16) 


(17) 


p q 

where  r is  the  distance  between  points  p and  Q and  n and  n are  the 

ri  p q 

outward  normals  to  the  surface  at  points  P and  Q,  respectively.  In  evaluat- 
ing K2  and  F2>  the  point  at  which  0^  = 0^  is  excluded  from  the  integration. 
Substituting  Eqs,  (15)-(17)  into  Eq.  (14)  gives 

L 


J *<V  {mv  +w}  d‘ 
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l 

-*(V|  {VW  + w}  % 

t <18> 

- J v(V  {Il(rN)  + X2  (rpq)}  dSq 

o 

= 2 rr  j $(s  ) + cyv(s  ) j 
L p P J 

where  t is  the  length  of  the  generating  line  of  the  surface  of  revolution. 

The  s-0  coordinate  directions  have  now  been  essentially  uncoupled  so  that 
the  problem  has  been  reduced  to  the  evaluation  of  the  line  integrals  in 
the  coordinate  directions  on  the  surface  of  the  body.  This  formulation 
does  not  restrict  the  form  or  type  of  boundary  conditions  on  the  body; 
it  merely  assumes  that  the  boundary  conditions  can  be  represented  by  a 
sum  (expanded  in  a set)  of  tangential  modes. 

III.  Results 

The  acoustic  fields  for  a sphere,  cylinder,  and  inlet  configuration 
have  been  computed  by  numerical  solution  of  Eq.  (18)  using  the  techniques 
described  in  Ref.  10.  Basically,  this  method  consists  of  first  specifying 
the  p-z  coordinates  and  the  normal  vector  at  each  point  on  the  surface. 

From  these  quantities  the  distances  and  the  normal  derivatives  can  be  ob- 
tained. Tire  integral  in  Eq.  (18)  is  then  separated  into  n integrals  taken 
over  subintervals  of  length  <t/n.  The  acoustic  potential  is  assumed  constant 
over  each  subinterval  and  the  integrations  are  performed  numerically  using 
Gauss -Legend re  quadrature  in  the  p-z  plane.  Over  the  subinterval  containing 
the  point  P,  the  integrand  in  Eq,  (18)  becomes  infinite  since  r approaches 
zero.  Thus,  only  an  even  number  of  points  is  used  in  the  quadrature  algorithim, 

since  an  odd  number  would  necessitate  inclusion  of  the  point  where  r _ = 0. 

pq 


* 
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J 


A Gauss-Legendre  quadrature  formula  is  used  in  the  circumferential  direction 
to  evaluate  Eqs.  (15)  - (17).  All  calculations  were  performed  on  the  Georgia 
Tech  CUC  Cyber  70/7|5  with  sixteen  significant  figures. 

In  all  geometries  investigated,  exact  solutions  were  obtained  for  m = 0 
by  assuming  a monopole  source  located  at  point  ( p ,z)  = (0,0)  inside  the  sur- 
face. The  normal  velocities  and/or  admittance  values  are  then  computed  at 
each  point  on  the  surface  using  Eq.  (6)  and  taken  as  the  boundary  conditions 
in  Eq.  (18).  The  surface  potential  $ (s^l  is  then  computed  from  Eq . (18)  and 
the  far  field  potential  is  obtained  by  numerically  solving  Eq . (2)  with  Eq . (6). 
The  computed  surface  and  far  field  potentials  are  then  compared  with  the  known 


potential  distribution  of  the  monopolc  source 

ika 

CP(P)  - “ 


(19) 


whore  a is  the  distance  from  the  source  to  the  observation  point.  For  m = 1 
a dipole  source  was  used  to  generate  exact  solutions,  and  for  m = 2 a quadra- 
pole  source  was  used. 

To  investigate  the  effect  of  the  coupling  constant  a in  Eqs.  (15)  - (17), 
the  surface  potential  distributions  for  a sphere  of  unit  radius  with  a uni- 
formly vibrating  surface  (i.e.  m = 0)  were  computed  for  or  =0 , i , and  i/k.  Twenty 
subintervals  were  taken  in  the  p -z  plar.e,  a four-point  Gauss-Legendre  quadra- 
ture formula  was  used  over  each  subinterval  and  a twenty-point  Gauss-Legendre 
formula  was  used  in  the  0 direction.  The  magnitude  of  the  potential  should  be 
unity  at  all  points  on  the  surface.  The  results  presented  in  Fig.  3 show  the 
computed  magnitudes  of  the  surface  acoustic  potential  to  be  in  error  by  12 
per  cent  for  a - 0 at  nondimensional  wave  numbers  ka  close  to  tt  i 2 n > and  3n. 

These  results  are  those  that  would  be  obtained  from  Eq.  (7).  The  relatively 

l. 

large  errors  are  expected  from  the  analysis  of  Burton  and  Miller  and  from  previous 

2,5 

investigations  using  Eq.  (7).  Burton  proves  that  setting  the  imaginary 

r 
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part  of  a nonzero  guarantees  unique  solutions  to  Eq.  (18),  For  a - i the 
maximum  error  is  reduced  to  less  than  four  per  cent  except  when  k is  close 
to  8.0.  However,  when  o=  i,  and  for  sufficiently  high  values  of  ka,  Eq.  (9) 
is  dominated  by  terms  arising  from  Eq.  (8).  As  a result,  the  solution  equa- 
tions become  ill-conditioned  when  ka  is  sufficiently  high  and  close  to  one 
of  the  eigenf requencies  associated  with  the  integral  equation  based  on  Eq . (8). 
In  'fable  I computed  values  close  to  these  eigenf requencies  and  the  eigenfre- 
quencies  of  Eq . (7)  are  compared  with  exact  results  for  q 0,  i,  and  i/k. 

In  all  cases,  the  value  of  i/k  gives  the  most  accurate  results.  In  Table  II, 
the  effect  of  introducing  an  admittance  condition  is  presented  for  or  = i/k. 

The  admittance  Y(Q)  and  forcing  function  A(Q)  in  Eq.  (6)  are  chosen  so  that 
the  relations 

fj  - *(Q)cp  = A(Q)  ; cp  = ~~~  (20) 

q 

are  satisfied  on  the  surface  and  the  exact  solutions  can  be  readily  computed. 
The  loss  in  accuracy  when  an  admittance  condition  is  used  is  minimal  and  re- 
stricted to  the  third  significant  figure.  However,  for  discontinuous  boundary 
conditions,  where  the  forcing  function  is  specified  over  one  part  of  the  sur- 
face (i.e.,  the  admittance  is  zero  there)  and  the  admittance  is  specified 
over  the  remaining  surface,  errors  of  over  ten  per  cent  in  the  real  and  imagi- 
nary parts  of  the  computed  surface  potential  result.  For  comparison,  the  case 
of  a constant  forcing  function  and  admittance  over  the  sphere  for  Of  = 0 is  also 
presented  and  in  all  cases  yields  results  of  less  accuracy  than  those  obtained 
with  a ~ i/k. 

Tn  this  study  consistently  good  results  were  obtained  with  a s i/k.  In 
Fig.  3 the  computed  and  exact  values  for  or  = i/k  agree  to  three  significant 
figures  over  the  range  of  nondimensional  wave  numbers  from  one  to  ten.  In 
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fact,  for  this  value  of  a , the  accuracy  is  significantly  better  at  all  wave 
. 4 

numbers  investigated.  While  Burton  and  Miller  provide  no  recommendations  for 

choosing  one  value  of  a over  any  other  value  with  an  imaginary  component, 

the  choice  o'  = i/k  used  in  the  present  study  can  be  explained  as  follows. 

2 

The  terms  in  Eqs.  (15)  - (17)  which  involve  a are  of  order  k whereas  the 

remaining  terms  are  of  order  k.  Therefore,  at  higher  wave  numbers  the  terms 
2 

of  order  k dominate.  By  choosing  a to  vary  inversely  with  the  wave  number, 

all  terms  in  Eqs.  (15)  - (17)  remain  of  the  same  order  with  respect  to  wave 
* 

number . 

A proo  iem  ol  iuo l c piTuctiCal  imp  ortance  is  the  finite  axisymmetric  duct 
since  this  surface  approximates  an  engine  configuration.  The  surface  poten- 
tial distributions  are  presented  in  Fig.  4 at  different  nondimensional  wave 
numbers  for  m = 0.  The  normal  acoustic  velocity  distribution  A(Q)  is  chosen 
so  that  the  solution  for  the  acoustic  potential  satisfies  Eq . (19).  The 
parameter  a is  taken  to  be  i/k.  Twenty  subintervals  are  taken  in  the  p -z 
plane  and  a twenty-point  Gauss-Legendre  quadrature  is  used  in  the  0 direction. 
Tn  Fig.  4 the  variations  of  the  magnitude  and  phase  with  distance  along  the 
perimeter  s are  presented.  The  largest  errors  in  the  computed  magnitude  of 
the  potential  of  about  ten  per  cent  occur  on  the  ends  of  the  cylinder  and  at 
the  corners.  The  results  at  the  ends  can  be  improved  without  increasing  the 
number  of  points  by  area  weighting  rather  than  taking  equidistant  points 
along  the  perimeter.  The  errors  in  the  phase  are  less  than  four  per  cent 
in  all  cases.  The  errors  in  the  magnitude  of  the  computed  surface  potential 
increase  with  increasing  nondimensional  wave  number; 

* It  is  interesting  to  note  that  in  the  report  by  Chertock^  he  suggests  the 
use  of  1/k  on  the  grounds  that  it  has  the  correct  physical  dimensions 
(i.e.  length)  that  will  maintain  the  dimensional  homogeneity  of  Eq.  (18). 
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but,  even  when  ka  = 10,  the  numerical  results  are  within  ten  per  cent  of 
the  exact  solutions.  For  0 = 0 or  i the  errors  are  significantly  larger 
above  ka  = 2. 

In  most  inlet  problems  the  boundary  conditions  are  discontinuous  with 
the  acoustic  velocity  or  potential  (which  is  directly  proportional  to  die 
acoustic  pressure)  specified  over  part  of  the  surface  and  the  admittance 
(representing  liners)  over  the  rest.  To  determine  the  effect  of  the  discon- 
tinuities and  the  use  of  an  admittance  function  on  the  numerical  results 
for  m = 0,  a cylinder  was  investigated.  The  velocity  was  specified  on  the 
ends  and  the  admittance  was  specified  in  the  center  so  that  the  solution 
for  § was  given  by  Eq.  (19)  and  Eq.  (6)  is  satisfied.  Again,  twenty  points 
are  used  in  the  p-z  and  Q directions.  The  results  are  shown  in  Fig.  5. 

Although  the  errors  in  the  numerical  results  for  this  case  are  higher  than 
those  observed  in  Fig.  4,  the  errors  still  remain  within  10  per  cent  for 
values  of  ka  less  than  5.  However,  when  ka = 10  errors  of  up  to  40  per  cent 
in  the  magnitude  of  the  potential  are  encountered  close  to  the  discontinuity 
in  the  boundary  condition.  This  error  can  be  reduced  by  increasing  the  num- 
ber of  subintervals  in  the  p-z  plane.  Doubling  the  number  of  subintervals 
halves  the  error.  When  both  the  normal  acoustic  velocity  and  the  admittance 
are  continuous  on  the  surface,  the  errors  are  of  the  same  order  of  magni- 
tude as  those  of  Fig.  4.  For  tangential  modes,  the  variation  in  the  circum- 
ferential direction  behaves  as  cos  m 0 where  m = 0,1,2.....  To  check  the 
numerical  integration  scheme  rn  the  circumferential  direction,  the  surface  acous- 
tic potential  was  computed  for  m = 1 and  m = 2 for  the  cylinder  shown  in  Fig.  4. 
The  results  are  presented  in  Fig.  6 for  ka = 2 with  the  normal  acoustic  velocity 
specified  everywhere  on  the  surface.  The  computed  and  exact  results  (i.e.  from  a 
dipole  and  qudrapole)  are  in  agreement  to  within  two  per  cent  for  both  m = l 
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and  m = 2 , 

It  has  been  shown'’  that  once  the  surface  potential  has  been  accurately 
computed,  the  far  field  can  be  determined  to  at  least  the  same  accuracy  as 
the  surface  potential.  This  result  is  confirmed  by  the  data  presented  in 
Fig.  7 for  the  cylinder  of  Fig.  4 with  the  velocity  specified  everywhere  on 
the  surface  with  ka  = 2 and  m = 0.  The  results  at  20  radii  from  the  surface 
are  in  agreement  with  exact  results  obtained  from  Eq.  (19)  to  within  one 
per  cent  even  though  the  surface  errors  at  some  points  are  above  two  per 
cent.  Data  in  Fig.  8 show  that  accurate  results  are  obtained  at  distances 
greater  than  one  integration  stepsize  from  the  surface.  At  closer  distances 
errors  from  the  numerical  evaluation  of  the  singularity  in  the  Green's 
function  defined  by  Eq.  (4)  leads  to  large  errors. 

The  studies  of  the  acoustic  fields  of  the  sphere  and  cylinder  served 
to  check  out  and  refine  the  numerical  procedures  and  programming  techniques. 
The  next  configuration  investigated  was  an  inlet  used  in  a study  by  NASA.^ 
This  inlet  is  shown  in  Fig.  9 and  was  chosen  because: 

(1)  unlike  most  inlets  used  in  research  studies,  it  does  not  have  a 
bell-mouth  shape  but  is  shaped  like  a typical  inlet  used  in  exist- 
ing aircraft;  and 

(2)  complete  details  on  generating  the  inlet  boundary  are  given  in  Ref.  11. 
For  this  inlet,  all  cases  were  investigated  with  o = i/k. 

As  seen  in  Fig.  10,  the  normal  velocity  distribution,  which  represents 
a forcing  function,  is  highly  discontinuous  and  provides  a severe  test  of 
the  numerical  techniques  employed.  The  numerical  and  exact  solutions  for  the 
surface  acoustic  potential  are  compared  in  Fig.  10  for  32  and  54  subintervals 
taken  along  the  perimeter  of  the  inlet  in  the  p-z  plane.  Because  of  the  errors 
in  approximating  the  lengths  of  each  subinterval,  the  exact  solutions  differ 
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slightly  as  the  distance  along  the  perimeter  s increases.  The  centerbody 
in  Fig.  9 extends  from  0 < s < 0.8,  the  fan  inlet  covers  0.8  < s < 1.4, 
the  interior  contour  extends  from  1.4  < s < 3.5,  the  exterior  from  3.5 
< s < 5.5,  and  the  circular  arc  lies  within  the  interval  5.5  < s < 7.45. 
Increasing  the  number  of  points  decreases  the  error  proportionately  as 
indicated  by  the  data  in  Fig.  10  at  a nondimensional  wave  number  ka  of 
unity  where  a is  the  radius  of  the  inlet  at  the  fan  entrance  section. 

The  absolute  average  error  in  the  results  decreases  from  10.2  per  cent 
for  32  subintervals  to  4.16  per  cent  for  53  subintervals.  The  computation 
Lime  increased  from  53  seconds  to  143  seconds,  respectively. 

As  shown  in  Fig.  11,  the  errors  increase  with  increasing  frequency. 
Like  the  cylinder,  the  maximum  error  in  the  acousitc  potential  for  the 
inlet  configuration  occurs  at  the  points  of  discontinuity.  The  average 
error  increases  from  4.16  per  cent  at  ka  = 1 to  15  per  cent  at  ka  = 10. 

For  the  data  in  Figs.  10  and  11,  the  acoustic  potential  is  assumed 
constant  in  the  tangential  plane.  The  results  for  a cos(ra  0 ) distribution 
are  presented  in  Fig.  12  at  ka  = 2.  These  results  show  the  insensitivity 
of  the  accuracy  of  the  computed  results  to  the  tangential  distribution 
for  m ~ 1,2.  The  exact  soltuions  were  again  generated  by  assuming  dipole 
and  quadrapole  sources  located  at  ( p ,z)  = (0,0). 

Basedonthe  above  results  our  numerical  and  programming  techniques  are 
capable  of  yielding  reliable,  results  for  arbitrary  geometries  and  boundary 
conditions.  At  higher  frequencies,  (ka  > 5)  it  appears  that  more  points 
must  be  taken  to  increase  the  accuracy  of  the  computed  results. 


IV.  Summary  and  Conclusions 


An  integral  solution  of  the  Helmholtz  equation  is  developed  for  use  in 
acoustic  radiation  problems.  Unlike  the  classical  formulation  which  can  lead 
to  integral  equations  that  do  not  have  unique  solutions  at  frequencies  cor- 
responding to  certain  internal  eigenfrequencies  of  the  region  enclosed  by 
the  surface  under  consideration,  the  formulation  used  in  this  study  is  valid 
at  all  frequencies.  Also,  unlike  most  current  methods  and  formulations  it 
is  straight  forward  to  implement  regardless  of  how  complicated  the  surface 
or  the  boundary  conditions  may  be.  Hie  surface  potentials  computed  numeri- 
cally for  a sphere  and  cylinder  using  20  subintervals  along  the  perimeter 
and  for  an  inlet  configuration  with  53  subintervals  are  accurate  to  within 
ten  per  cent  for  nondimens ional  wave  numbers  ka  of  from  one  to  ten  where  k 
is  the  wave  number  and  a is  the  characteristic  length.  For  discontinuous 
boundary  conditions,  the  numerical  and  exact  values  are  in  agreement  to  with- 
in 10  per  cent  for  ka  < 5.  At  higher  frequencies  the  results  are  as  much  as 
40  per  cent  in  error  at  points  of  discontinuity  which  suggests  taking  more 
points  in  evaluating  the  integral  equation  to  increase  the  accuracy  when  dis- 
continuous boundary  conditions  are  specified.  Increasing  the  number  of  sub- 
intervals decreases  the  error  proportionately.  At  distances  greater  than  the 
numerical  integration  stepsize,  the  far  field  results  are  at  least  as  accu- 
rate as  the  corresponding  surface  potential  solutions. 
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Tabic  I 

Effect  of  the  coupling  parameter  a on  the  computed  values  of  the  surface 
potential  for  a sphere.  On  the  surface  A(Q)  = (l-ikje^  , Y(Q)  = 0,  <p 

CiXclC  t_ 

constant,  m = 0.  All  values  of  ka  correspond  to  internal  eigenfrequencies. 
Twenty  subintervals  were  taken  in  the  p-z  plane. 


EXACT 


0.998 

-0.012 


4.493409 


7.725252 


1 

0 

i/k 

2.0 

1.000 

-6.3 

0.001 

0.190 

0.217 

0.979 

0.976 

-2.0 

-1.000 

12.6 

0.000 

-0.081 

-0.128 

-0.994 

-0.992 

2.0 

1.000 

-19.0 

0.000 

0.308 

0.955 


-0.996 

0.031 


-0.400 

-0.872 


0.995 

-0.050 


0.217 

0.976 


-0.128 

-0.992 
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Table  II 

Effect  of  specifying  an  admittance  on  the  computed  surface  potential 


for  a sphere.  In  all  cases  m = 0,  twenty  subintervals  are  taken  in  the  p-z 

ilc  iic 

plane,  and  cPexl.t=*c  everywhere  on  the  surface.  For  Case  I,  A(Q)  = e 

j ilc 

and  Y(Q)  = .'A  everywhere  on  the  surface.  For  case  III,  A(Q)  = e (1-ik)  and 

Y(Q)  = 0 over  1/5  of  the  surface  and  A(Q)  = 0,  Y =-(l-ik)  over  the  remainder. 

Case  II  is  considered  in  Table  I. 


Effect  of  the  Coupling  Constant  on  the  Computed 
Surface  Potential  for  a Sphere  of  Unit  Radius  with 
20  Subintervals. 


Computed  Surface  Potential  for  a Cylinde 
First  and  Second  Tangential  Modes  for  ka 


R(F?Q)  =20 


Figure  7.  Computed  Far  Field  Potential  Distribution  for 
a Cylinder  at  k = 2,  m = 0,  and  20  Radii  from 
the  Center. 


o Computed,  m=1 


Appendix  C 


Paper  presented  at  the  AIAA  5th  Aeroacoust ics  Conference  in 
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Abs tract  exterior  sound  fields  to  be  a single  entity,  thus 

elininating  the  errors  associated  with  many  of  the 

Results  are  obtained  by  numerical  integration  previously  used  approaches.  Comparison  of  the  solu- 

of  a cy 1 indrical ly  symmetric  integral  represents-  tions  obtained  using  this  method  with  exact  solu- 
tion of  the  exterior  solutions  of  the  Helmholtz  Cions^  and  with  those  of  related  investigations 

equation  which  is  valid  (yields  unique  solutions)  should  shed  some  light  upon  the  applicability  of 

at  all  wavenumbers.  The  admittance  values  across  the  analytical  approaches  utilized  in  these  other 


the  entrance  plane  of  hard  walled  ducts  of  various 
lengths  and  geometries  are  computed  and  compared 
with  certain  "classical"  (e.g.,  Weiner-Hopf  ) va  lues  . 
The  internal  wave  structure  is  also  investigated 
for  straight  hard  walled  ducts  and  compared  with 
results  obtained  from  other  theories.  The  radiated 
sound  fields  from  ducts  of  different  geometries  are 
then  compared  for  both  unlined  and  Lined  configura- 
tions. It  is  found  that  changes  in  the  duct  geometry 
result  in  significant  changes  in  the  radiated  sound 
field.  Thus,  it  is  concluded  that  the  sound  suppres- 
sion by  liners  predicted  from  the  study  of  straight 
ducts  may  not  be  applicable  to  aore  complicated 
geometries  such  as  inlet  conf igurations. 

Int  reduction 

The  development  of  an  analytical  method  for  pre- 
dicting the  sound  field  radiated  from  axisymaetr ic , 
finite  length  ducts  is  of  much  practical  interest 
in  the  area  of  aeroacous tics , especially  for  the 
determiriat ion  of  the  sound  radiated  from  a turbofan 
inlet,  as  having  such  a capability  can  eliminate 
most  of  the  costly  full  scale  testing  presently  re- 
quired. In  a majority  of  past  investigations  of  the 
sound  radiated  from  ducts,  either  the  radiation 
problem  has  been  completely  ignored! or  the  duct 
acoustics  and  the  radiation  problem  have  been  treat- 
ed separately . 5 * In  these  studies,  the  behavior 
of  the  waves  inside  the  duct  was  determined  by  spec- 
ifying some  heuristic  boundary  condition  (e.g.,  a 
reflection  coefficient)  at  the  duct  entrance.  In  the 
latter  references  this  solution  was  then  used  to 
determine  the  sound  distribution  at  the  entrance 
plane  of  the  duct  and  this  was  used  to  predict  the 
properties  of  the  radiated  sound  field.  In  reality 
the  sound  fields  inside  and  outside  the  duct  are 
not  separate  entities  (i.e.,  they  are  coupled)  and, 
therefore,  they  cannot  be  properly  treated  separate- 
ly. These  "separate"  treatments  of  the  interior  and 
exterior  sound  fields  in  the  duct  radiation  problem 
undoubtedly  Introduce  errors  whose  determination 
requires  comparison  with  available  exact  solutions. 

In  this  paper  the  duct  sound  radiation  problem 
is  investigated  by  utiliiing  an  integral  solution? 
that  considers  the  combination  of  the  interior  and 
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investigations. 

The  solution  approach  utilized  in  this  study 
consists  of  the  numerical  solution  of  a special  in- 
tegral representation  of  the  solutions  of  the  Helm- 
holtz equation  for  an  exterior  (i.e.,  to  a given 
body)  domain.  The  applicability  and  accuracy  of  this 
solution  aporoach  have  been  demonstrated  for  two 
dimensions -0»! * , three  dimensional**, 12#  and  axisytu- 
netric5,13  geometries  in  earlier  investigations  con- 
ducted under  this  program  where  it  has  been  shown 
that  the  developed  approach  yields  unique  solutions 
at  all  wavenumbers. 

Solution  Procedure 

The  basis  of  this  method  is  set  forth  in  great 
detail  in  Ref.  14  and  therefore  will  not  be  repeated 
here.  In  related  studies  conducted  under  this  APOSR 
program  the  applicability  of  this  approach  to  the 
solution  of  acoustics  problems  involving  two  dimen- 
sional, three  dimensional,  and  axisymmetric  geome- 
tries has  been  demonstrated.  Since  the  analytical 
developments  and  resu2cs  of  these  studies  have  been 
published  elsewhere? »** , they  also  will  not  be  re- 
peated here.  Instead,  soae  of  the  advantages  of 
this  solution  approach  will  be  presented.  First,  as 
stated  earlier,  this  method  treats  the  duct  radia- 
tion problem  consisting  of  the  sound  generation, 
sound  propagation  and  reflection,  and  the  sound 
radiation  to  the  outside  as  a whole  without  sepa- 
rating it  into  its  component  parts  as  has  been  done 
in  related  irrves  ti  gat  ions.  Second,  the  method  can 
readily  handle  the  infinite  domains  encountered  in 
radiation  problems.  This  is  accomplished  by  employ- 
ing a fundamental  solution  G(P,Q)  which  satisfies 
the  Sooner field  radiation  conditions  in  the  inte- 
gral equation.  In  the  present  study  the  free  space 
Green's  function  has  been  used: 


ikr(P,Q) 
G(P*Q)  ' r(P,Q)  " 


where,  as  shown  In  Fig.  1,  Q is  a point  on  the  sur- 
face of  the  body  S,  P is  a point  in  the  exterior 
domain,  r(P,Q)  is  the  distance  between  these  points, 
and  k is  the  wavenumber.  The  third  advantage  is 
that  the  computer  program  developed  in  a related 
investigation9.13  is  quite  general  and  is  applic- 
able to  a variety  of  acoustic  radiation  problems 
involving  arbitrary  geometries  and  variable  bound- 
ary conditions.  This  computer  program  can  be  ap- 
plied to  different  problems  by  merely  changing  the 
input  data. 
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Fig.  1.  General  Geometry  of  the 
Radiation  ProbLem. 


It  has  been  shown  in  Ref.  8 that  unique  solu- 
tions of  the  external  acoustic  radiation  problem 
can  be  obtained  at  all  wavenumbers  by  solving  the 
following  integral  equation: 


It  is  shown  in  Ref.  9 that  maximum  accuracy  is  ob- 
tained when  cr  * i/k.  In  its  present  form,  Eq,  (2) 
contains  no  non-integrablc  singularities,  a fact 
which  considerably  simplifies  its  numerical  solu- 
tion. Thus  the  solution  of  the  entire  duct  radia- 
tion problem  has  been  reduced,  with  the  present 
formulation,  to  that  of  the  solution  of  a surface 
integral  over  the  body  S. 

The  problem  is  further  simplified  in  the  present 
study  by  Uniting  attention  to  axisymmetric  config- 
urations. In  this  case  Eq.  (2)  is  further  reduced 
coW 

j 4(Q)  (mm)  + K2(P,Q)}  dsq 
0 l 

- i(P)  J {Fj(P,Q)  + F2(P,Q)}dsq 

l ° (3) 

- V(Q)  {lx(P,Q)  + I2(P,Q)}  dsq 

- 2*{4(P)  + cr  V(P)} 


J f(«p(Q)  - G(P,Q)  dS 

V 3nq  anq 

+ a J J (w(Q)  - P(P))  L'^1^  dS 
V JV",  q 

(2) 

- a <s(P)  I f (n  • n )(ik)2  G(P,Q)  dS 
J J p q q 


where  s is  the  distance  _ long  the  2-D  projection  of 
the  body  in  the  p-z  plane  (i.e,.  See  Fig.  2),  and 
the  influence  functions  and  I2  are  given  by 
TT 

^(P.Q)  « 2 G(P,Q)  (coa  m eq)  d0q 

n (4) 

I2(?.Q)  - 2a  jo  (cos  m Gq>  d0q 


- 

Snp  1 


2 tt(<p(P)+ ) 


where,  as  shown  in  Fig.  2,  the  point  P has  been 
moved  to  the  surface  of  the  body  S,n  represents  an 
external  normal  from  the  body,  represents  a 

normal  derivative  (vn),  and  a is  a complex  coupl ing 
cons  tant. 


the  kernel  functions  and  K2  are 
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Fig.  2.  General  Aaclf yuroetrlc  Geometry 
and  2-0  Projection. 


the  forcing  functions  Fj  and  F2  are 

FjCP.Q)  - 2o  ^ C(P,Q)  (ik)2(np«nq)d0q 


and 


(6) 


Vp.<»  ■ 2“  £ d0q  \*% 


- V 
on 


cos  m 0 




(7) 


cos  ra  0 


following  the  notation  of  Fef.  15.  In  the  above 
notation  m is  the  tangential  mode  and  0^  has  been 
taken  as  zero  (i.e.,  cos  m 0 « l). 


Introducing  ax isyranetry  further  simplifies  the 
solution  of  Eq.  (2)  to  the  evaluation  of  line  in- 
tegrals on  the  2-D  projections  of  the  body  S in  the 
p -Z  plane  as  shown  in  Fig.  2.  Also,  the  formula- 
tion is  valid  for  all  tangential  modes;  however, 
each  mode  must  be  solved  for  separately. 


The  required  soLution  is  obtained  by  solving 
Eq«  (3)  for  the  surface  distribution  of  the  acoustic 
potential  $ or  the  normal  acoustic  velocity  V, 
whichever  is  unknown.  Also,  solutions  can  be  ob- 
tained for  $ by  using  an  effective  admittance,  de- 
fined as  Y -V/$  » do/dn/p,  as  a boundary  condition 
over  any  part  of  the  surface  of  the  body  (e.g., 
replace  V by  4>Y  in  Eq.  (3)  at  the  points  where  Y is 
known  on  the  bodyt ) # 

Once  the  surface  distributions  of  the  acoustic 
potential  and  the  normal  acoustic  velocity  are 
known  on  the  surface  of  the  body,  the  radiated 
sound  field  can  be  determined  using  the  following 
integral  representations1^  for  cp  and  ^ in  the 
field:  dn 

SJ  - c(p>Q>  *fP>dv  4r*(p> 

(8) 

and 

r r,_(Q)  afciLfti . mlsi  ^m)ds , ^ Mil 

JSJ  SVn4  a"p  8nq  q anp 

(9) 

where  the  point  P is  now  located  in  the  space  sur- 
rounding the  body  (i.e.,  See  Fig.  1 ) so  that  the 
kernel  functions  are  no  longer  singular.  For  axisy- 
mmetric  bodies,  Eqs.  (8)  and  (9)  reduce  to: 

J ($(Q)  K,(P,Q)  - In  (P,Q)  V (Q) )ds  - 4n  *(P) 

q (10) 

and 

f*' 

J <*<<?)  K2(PpQ)  - V (Q))ds  - 4tt  V(P) 

4 (11) 


Results 


In  the  present  investigation,  the  integral 
solution  technique  has  been  utilized  to  study  the 
dependence  of  the  radiated  sound  field  and  the 
acoustic  characteristics  of  the  duct  upon  the  duct 
geometry  and  the  acoustic  properties  of  the  duct 
wall.  This  investigation  has  been  carried  out  with 
the  objective  of  evaluating  the  dependence  of  the 
sound  field  radiated  from  a jet  engine  inlet  on 
the  inlet  characteristics,  and  to  evaluate  the 
validity  of  the  analytical  approaches  and  assump- 
tions utilized  in  related  investigations.  First, 
the  effect  of  the  length  of  the  duct  on  the  ad- 
mittance values  at  the  entrance  plane  of  the  duct 
was  investigated.  The  purpose  of  this  study  was 
not  only  to  show  the  effect  of  changing  the  L/a 
(i.c.,  length/radius)  of  the  duct  but  also  to  show 
that  the  admittance  is  not  a constant  across  the 
entrance  plane.  This  is  significant  as  many  in- 
vestigations of  similar  and  related  problems  assume 
the  existence  of  a constant  "reflection  coefficient" 
at  this  plane.  Therefore,  the  effect  of  the  radi- 
ated sound  field  on  the  duct  acoustics  cannot  be 
properly  accounted  for.  Second,  the  effect  of  the 
internal  geometry  of  the  duct  on  the  radiated 
sound  field  was  investigated.  This  was  done  to 
show  that  changes  In  the  Internal  geometry  of  the 
duct  result  in  large  changes  in  the  admittance 
values  at  the  duct  entrance  and,  therefore,  In 


large  changes  in  the  radiated  sound  field;  a fact 
which  is  often  ignored  in  related  studies.  Third, 
the  internal  wave  structure  in  the  duct  was  investi- 
gated to  show  that  th  is  integral  solution  technique 
can  predict  both  the  internal  and  external  sound 
fields.  In  this  connection  it  should  be  pointed  out 
again  that  this  solution  technique  automatically 
accounts  for  the  coupling  effects  between  the  in- 
ternal and  external  sound  fields.  Fourth,  the  radi- 
ated sound  fields  for  two,  dimensionally  similar, 
acoustically  lined  duct  configurations  (i.e.,  a 
straight  duct  and  an  engine  inlet)  are  compared  to 
show  that  optimum  admittance  values  for  liners,  de- 
termined from  the  study  of  their  effectiveness  in  a 
straight  duct,  do  not  necessarilly  carry  over  to  the 
core  complicated  inlet  configurations. 

In  this  study,  the  surface  distributions  of  the 
unknowns  of  the  problem  (i.e.,  the  acoustic  poten- 
tial and/or  the  normal  acoustic  velocity)  are  ob- 
tained from  the  numerical  solution  of  Eq.  (3).  Since 
the  unknown  functions  <?(Q)  or  V(Q)  appear  in  the  in- 
tegrands, it  is  necessary  to  solve  a square  matrix; 
thus,  the  re  uired  computing  time  increases  roughly 
as  the  square  of  the  number  of  points  taken  on  the 
surface  of  the  body  S.  In  a typical  run  on  the 
Georgia  Tech  CDC  Cyber  70/74  computer, 140  seconds  of 
computing  time  were  required  to  solve  for  the  sur- 
face values  of  either  ♦ or  V at  53  points  on  the 
body. 

To  obtain  values  for  $ and  V in  the  far  field, 
Eqs.  (10)  and  (11)  are  solved  by  simple  numerical 
integration.  The  time  required  for  this  computation 
is  roughly  proportional  to  the  number  of  points  in 
the  field  and  the  number  of  points  on  the  body.  In 
a typical  run,  70  seconds  of  computing  time  was  re- 
quired to  calculate  both  $ and  V at  57  points  in 
the  field  with  53  points  on  the  body. 

To  determine  the  dependence  of  the  admittance 
at  the  duct  entrance  plane  on  the  duct  length,  these 
admittance  values  were  computed  for  hard  walled 
straight  ducts  having  different  L/a  values,  as  shown 
in  Fig.  3. 


Fig.  3.  Geometry  Used  for  Straight  Duct 
Computations. 


In  this  case,  driving  consisted  of  a unit  acoustic 
velocity  (V*»l)  across  the  driver  face.  The  com  - 
puted  admittance  values  at  the  duct  entrance  plane 
are  plotted  for  two  non-dimensional  wavenumbers 
(i.e.,  ka  - 1 and  3)  in  Figure  4.  Also  noted  for 
comparison  are  the  "classical"  values  for  flanged 
and  unflanged*^  pipes.  These  results  indicate  that 
even  for  plane  wave  sound  excitation  the  admittance 
at  the  duct  entrance  plane  varies  with  the  trans- 
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Pig.  4.  Admittance  at  the  Exit  Plane  of  a 
Straight  Duct  vlth  the  Classical 
Values  for  Flanged  anJ  Unflanged  Pipes. 


To  determine  geometrical  effects,  the  admittance 
at  the  entrance  plane  of  a hard  vailed  Inlet  config- 
uration^® (l.e.,  See  Fig.  5.),  vlth  L/a  - 2 was  also 
calculated  assuming  the  same  type  of  source  excita- 
tion for  comparison.  These  results  are  plotted  In 
Figure  6.  When  compared  with  the  corresponding  re- 
sults for  a straight  duct  of  the  same  basic  dimen- 
sions (See  Fig.  4.)  it  is  seen  that  the  admittance 
values  change  significantly.  This  Is  but  one  Indica- 
tion of  the  Importance  of  the  need  to  properly 
account  for  the  Internal  geometry  of  the  duct  vhen 
investigating  duct  radiation  problems. 

The  Internal  wave  structure  was  also  Investi- 
gated for  two  hard  vailed,  straight  ducts  to  de- 
termine its  dependence  upon  the  characteristics  of 
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Fig.  6.  Admittance  at  the  Exit  Plane  of 
the  Inlet  (L/a  - 2) 


4 


the  sound  source.  The  resul's  ate  presented  in 
Fig.  for  two  different  drivers,  one  having  a 
constant  unit  normal  acouscic  velocity  (i.e.,  V-l) 
and  the  other  a cosine  distribution  (i.e.,  V=  cos 
(rtp)  where  OSpSl),  at  the  sound  source  plane, 
for  the  case  where  ka  = land  m (the  tangential  mode 
number)  equals  zero.  These  results  show  that  the 
difference  between  Che  sound  sources  quickly  dis- 
appear with  increasing  distance  from  the  source 
plane,  as  expected  for  low  values  of  the  non-dimen- 
sional wavenumber  ka.  The  results  for  two  addition- 
al cases  are  also  presented; 
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Fig.  8.  Radial  Distribution  of  the  Acoustic 
Potential  in  a Hard  Ualled  Straight 
Duct  (L/a  = 3,  ka  «=  3,  m » 0) 
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Fig.  7.  Radial  Distributions  of  the  Acoustic 
Potential  In  a Hard  Walled  Straight 
Duct  (L/a  ■ 3,  ka  • 1 , m ■ 0) 


in  Fig.  8 the  same  two  sound  sound  sources  are  used 
for  the  higher  wavenumber  of  ka  - 2 while  in  Fig.  9 
a sound  source  with  a normal  acoustic  velocity  dis- 
tribution of  V - sin  (ptr)  and  a tangential  mode 
number  m » 1 Is  used.  All  of  these  results  indicate. 
In  agreement  with  basic  acoustics,  that  the  finer 
details  of  the  sound  source  are  "washed  out”  within 
a length  on  the  order  of  ka  from  the  sound  source 
plane , 


Finally,  the  radiated  sound  fields  from  both  a 
straight  duct  and  an  inlet  configuration  were  com- 
puted and  compared  for  both  lined  and  unlined  walls. 
The  values  used  for  the  wall  admittances  and  the 
wavenumbers  were  chosen  to  be  the  same  as  thos< 
used  by  Zorumski  in  Ref.  19  so  that  comparisons 
could  be  made.  Two  different  liners  were  run;  a 
constant  admittance  liner  and  a segmented  liner. 

In  the  case  of  Che  straight  duct  the  results  do  not 
show  Che  same  difference  in  the  radiated  sound  pres- 
sure level  between  the  segmented  and  constant  ad- 
mittance liners  (See  Fig.  10.)  as  did  Zorumski  (See 

Fig.  11.). 
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Fig.  9.  Radial  Distribution  of  the  Acoustic 
Potential  in  a Hard  Walled  Straight 
Duct  (L/a  = 3,  ka«l,m  = l) 
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Fig.  11.  Far  Field  Radiation  From  A Straight 
Duct  as  Per  Ref*  19  (L/a  « 2). 


To  investigate  the  effect  of  geometry,  the  inlet 
configuration  shown  in  Fig.  5 was  run  with  the  same 
two  liners  at  the  same  non-dimensional  wavenumbers. 
The  opposite  trends  to  those  predicted  by  Zorumski 
were  noted  (See  Fig.  12.);  that  is,  that  the  seg- 
mented liner  is  superior  to  the  constant  admittance 
liner  in  reducing  the  radiated  sound  pressure  levels. 
Thus,  it  is  concluded  that  geometrical  details  can 
significantly  affect  the  characteristics  of  the 
sound  power  radiated  from  dimensionally  similar  ducts 
(i.e,,  ducts  having  the  same  L/a  ratios.). 
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Determination  of  the  Adnittance  of  the  Liner 
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Nomenclature 

A orifice  area 

c speed  of  sound 

d orifice  diameter 

f frequency 

fQ  resonant  frequency  of  the  resonator 

k wave  number 

L backing  depth 

^EFF  effective  orifice  length 

n normal 

p acoustic  pressure 

t orifice  length 

V cavity  volume 

v acoustic  velocity 

Y specific  admittance  of  liner 

y effective  admittance  of  liner 

Z specific  impedance  of  liner 

a absorption  coefficient 

o open  area  ratio  of  the  liner 

o open  area  ratio  of  a resonator 

p density 

ll  coefficient  of  viscosity 

tt  3.1415926 

cp  acoustic  potential 

0 specific  resistance  ratio 

X specific  reactance  ratio 


40 


The  liner,  which  was  employed  in  the  preliminary  testing  for  this 
AFOSR  contract,  was  originally  developed  for  another  program.  Since  it 
consists  of  a matix  of  Helmholtz  resonators  it  is  highly  tuned;  that  is,  it  has 
a rather  sharp  absorption  peak.  As  discussed  in  Section  II,  the  liner  was 
designed  for  use  above  the  IT  cut-off  frequency  and  thus  shows  little 
absorption  below  this  frequency  where  most  of  the  runs  for  this  program  will 
be  made.  Thus  the  liner  must  be  retuned  so  that  it  becomes  effective  in  the 
frequency  range  where  most  of  the  testing  will  be  done.  In  this  appendix 
basic  liner  theory  is  reviewed  and  a relatively  simple  redesign  of  the  liner  is 
proposed  to  make  it  more  effective  in  the  frequency  range  of  interest  in  this 
research  program. 

For  input  into  the  computer  program  certain  values  must  be  known. 
First,  the  sound  pressure  level  at  the  "driver  plane"  must  be  known. 
Experimentally  this  corresponds  to  the  nozzle-liner  plane  and  can  be  directly 
measured.  The  other  value  which  must  be  known  is  the  effective  admittance 
of  the  liner  defined  as 


y 


§2.  / 
Sn  / ?p 


(D-l) 


where  cp  is  the  acoustic  potential  and  ^ is  the  normal  acoustic  velocity 
defined  with  an  outward  facing  normal. 

In  Reference  I equations  are  given  to  calculate  the  specific  acoustic 
impedance  of  an  array  of  helmholtz  resonators,  that  is 

Z = 0 - i x (D-2) 


Li, 
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where  the  specific  resistance  ratio  0 and  the  specific  reactance  ratio  x 
are  defined  as 


(tW)172  (1  + l/d) 

ape 


2 nfo^EFE  f 
So  'fo 


For  the  definitions  of  the  variables  used  here  see  Figure  D-l.  It  will  be  noted 
here  that  these  definitions  assume  an  inward  facing  normal.  The  resonant 
frequency  of  the  resonator  is  defined  as 


o 2 n V V'C'Epp  2 n V Lc 


and  the  effective  orifice  length  is  found  to  be 


= t + o.85  d ( i - o.?yT ) 


The  specific  acoustic  admittance  is  defined  as 


* = \ - V/P 


and  since  the  acoustic  potential  and  the  acoustic  pressue  are  related  by 


p = i p c k cp 


the  effective  admittance  is  found  to  be 


y = -ipc  k Y = -ipc  k -Y'V  (D-8) 

(0  +X.  ) 

where  the  minus  sign  is  due  to  the  switch  from  inward  to  outward  facing 
normal  for  the  acoustic  velocity. 

The  liner  is  most  effective  in  damping  acoustic  waves  at  its  resonant 
frequency  which  for  our  case  is  ~ 740  Hz.  This  is  above  the  IT  mode  of  the 
duct  which  is  ~ 695  Hz  and  thus  the  full  effectiveness  of  this  liner  is  never 
achieved.  A plot  of  the  absorption  coefficient  of  the  liner,  a , vs.  frequency 
is  presented  for  the  present  liner  configuration  in  Figure  D-2.  Here  the 
absorption  coefficient  is  defined  as 

“ = — T (D-9) 

(0  + 1)  + x 

There  are  two  paths  that  may  be  taken  to  obtain  test  results  at 
maximum  liner  efficiency.  The  first  method  is  to  purposely  drive  a IT  wave 

O 

by  using  two  drivers  driving  180  out  of  phase.  This  will  create  a transverse 
acoustic  wave  structure  in  the  tube  and  the  results  can  then  be  compared  by 
using  a mode  number  of  1 (instead  of  zero)  in  the  computer  programs.  The 
other  alternative  is  to  redesign  the  liner  such  that  its  resonant  frequency  is 
reduced  below  the  IT  mode  cut-off  frequency.  This  can  be  accomplished  by 
increasing  the  backing  distance,  L.  By  increasing  L from  0.5"  to  0.775"  the 
resonant  frequency  of  the  liner  drops  to  ~ 685  Hz,  below  the  IT  mode.  If  the 
diameter  of  the  backing  cavity  is  also  increased  from  1.0"  to  1.25"  the 


resonant  frequency  of  the  liner  can  be  further  reduced  to~55 0 Hz.  Thus,  by 
simply  drilling  out  the  backing  cavity  the  tuning  frequency  of  the  liner  can 
be  altered  enough  so  that  its  absorption  peak  is  well  below  the  IT  cut-off 
frequency;  that  is,  in  the  range  of  frequencies  that  will  normally  be  used  for 
testing  in  this  research  program  (i.e.,  300-650  Hz). 
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When  considering  the  radiation  from  an  open  duct  it  is  found  that 
some  of  ..he  energy  is  radiated  and  some  reflected  (with  a phase  shift)  back 
down  the  duct.  It  is  common  to  associate  an  effective  admittance 
(impedance  or  reflection  coefficient)  at  the  exit  plane  with  this 
phenomonon.  In  this  note  three  methods  for  obtaining  an  effective 
admittance  are  compared. 

The  first  method  follows  the  analyses  of  Helmholtz  and  Rayleigh  in 
which  the  end  of  the  duct  is  approximated  by  a piston  radiating  into  a half 
space  from  an  infinite  baffle.  In  this  analysis  the  classical  integral 
representation  of  the  solutions  of  the  Helmholtz  equation  is  solved  with 
certain  approximations.  Results  for  this  configuration  (commonly  known  as  a 
flanged  pipe)  using  this  method  are  presented  in  Reference  1. 

The  second  method  consists  of  the  solution  of  a Weiner-Hopf  type 
integral  equation.  In  this  analysis  the  duct  is  assumed  to  be  semi-infinite  in 
length  and  infinitely  thin.  Results  for  this  type  of  analysis  are  presented  in 
Reference  2.  This  configuration  is  commonly  known  as  the  unflanged  pipe. 

It  is  interesting  to  note  that  these  two  configurations  represent  the 
logical  limits  of  this  type  of  problem  in  that  the  first  can  be  viewed  as  an 
infinitely  thick  duct  while  the  second  is  infinitely  thin.  Neither  of  these 
configurations  can  account  for  the  case  of  a duct  of  finite  length,  however. 

The  third  method  employs  a special  cylindrically  symmetric  integral 
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representation  of  the  exterior  solutions  of  the  Helmholtz  equation.  ’ ' 
Using  this  method  it  is  possible  to  calculate  the  acoustic  pressure  and 
velocity  anywhere  in  the  external  field  - including  the  inside  of  the  duct 
itself.  From  these  values  inside  the  duct  an  effective  admittance  can  be 
calculated  using  a simple  standing  wave  analysis  like  that  used  in  an 
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impedance  tube.  All  that  is  required  for  this  method  is  knowledge  of  two 
complex  acoustic  quantities  (i.e.,  amplitude  and  phase)  at  two  points  in  the 
duct. 

The  geometric  restrictions  on  the  integral  technique  are  that  the  duct 
(or  any  radiating  body)  must  be  a finite  (i.e.  no  infinitely  thin  walls)  closed 
body.  A sketch  of  the  duct  used  in  this  analysis  is  presented  in  Figure  1. 

To  determine  if  the  length  of  the  duct  makes  a significant  difference 
in  the  admittance  values  at  the  duct  entrance  computer  analyses  of  ducts  of 
varying  lengths  were  run  (i.e.  L/a  = 1,2,3)  at  two  different  non-dimensional 
wave  numbers  (i.e.  ka  = 1 and  3).  Since  the  method  is  capable  of  calculating 
the  actual  radially  varying  admittance  across  the  exit  plane  of  the  duct, 
these  are  presented  in  Figure  2.  The  driver  consisted  of  specifying  a unit 
normal  acoustic  velocity  while  on  the  rest  of  the  body  the  admittance  was 
specified  as  zero.  As  can  be  seen  the  length  of  the  duct  L has  little  effect  on 
the  admittance,  defined  as  the  ratio  of  the  component  of  the  acoustic 
velocity  normal  to  the  surface  to  the  acoustic  pressure.  Also  noted  for  the 
sake  of  comparison  are  the  values  for  the  flanged  and  unflanged  pipe  at  the 
appropriate  values  of  ka. 

Computer  analyses  of  a duct  with  L/a  = 3 were  then  done  with  the 
same  boundary  conditions  specified  as  above  for  various  values  of  ka.  For 
each  case  the  acoustic  potential  and  velocity  were  calculated  at  1 1 equally 
spaced  points  along  the  centerline  of  the  duct  from  Z = 1 to  Z = 2.  A 
standing  wave  analysis  was  then  done  employing  a Least-Square  method  to 
solve  the  overdetermined  system  of  equations.  The  results  of  these  analyses 
are  presented  in  Figure  3 along  with  the  values  for  a flanged  and  an 
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unflanged  pipe.  Similar  analyses  were  done  off  the  centerline  of  the  duct; 
however,  no  significant  differences  in  the  computed  values  for  the 
admittance  at  the  entrance  plane  of  the  duct  were  found.  Also  noted  for 
comparison  in  Figure  3 are  the  values  of  the  admittance  calculated  on  the 
centerline  of  the  duct  at  the  exit  plane. 
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